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FURTHER  STUDIES  ON 
NONLINEAR  ADAPTIVE  OPTICS 


by: 


A.  Elci,  J.  Nagel  and  D.  Rogovin 


Science  Applications,  Inc. 

1200  Prospect  Street 
La  Jolla,  California  92038 


Abstract:  This  work  examines  interband  transitions  near  the  bandedge,  alkali  type 
ious  embedded  in  dielectrics,  and  two-photon  resonances  in  three-level  systems 
for  the  purpose  of  four-wave  mixing  and  phase  conjugation  at  the  iodine  laser 
frequency.  It  also  examines  the  time-dependent  properties  of  unstable  resonators 
with  a  saturable  gain  medium.  The  third  order  susceptibility  and  the  conjugate 
reflectivity  of  semiconductors  are  calculated.  Near  the  band  edge  the 
reflectivity  is  enhanced.  The  unstable-resonator  field  is  represented  in  terms 
of  a  series  of  empty  resonator  modes,  and  a  set  of  first  order  in  time 
differential  equations  for  the  exapansion  coefficients  is  achieved.  There  is  a 
stable  solution  and  a  bistable  one,  with  self-sustaining  pulsations.^ 


1 


I.  STATEMENT  OF  WORK 


a.  Examine  non-linear  materials,  particularly  liquids  and  gasses  for 
candidates  for  high  efficiency,  degenerate,  four-wave  mixing  media  for  use  with 
an  Iodine  (I*)  laser.  Several  such  candidates  are: 

1.  Nd-YAG 

2.  I* 

3.  Two-photon  resonant  dyes 

4.  Two-Photon  resonant  metal  vapor 

5.  Excited  state  alkalai  metal 

b.  Consider  the  above  materials  and  such  others  that  may  arise  In  tne 
study,  or  be  suggested  by  the  Air  Force  Weapons  Laboratory. 

c.  Consider  systems  that  might  emulate  a  potential  non-linear  medium 
at  a  frequency  that  1s  experimentally  attainable,  I.e.,  a  D.F.  frequency. 

d.  Modeling  of  the  media  selected  will  be  developed  in  conjunction 
with  experiments  at  the  Air  Force  Weapons  Laboratory. 

e.  Analyze  the  effects  of  transverse  mode  competition  In  unstable 
resonator  lasers  using  Lamb  Theory  of  the  Active  medium  and  actual  unstable 
resonator  mode  functions. 


II.  NONLINEAR  ADAPTIVE  OPTICS  SUMMARY 


The  results  of  our  work  on  nonlinear  adaptive  optics  are  given  in  two 
papers:  "Novel  Approaches  to  Nonlinear  Adaptive  Optics  at  1.3u",  which  will  be 
published  in  the  Proceedings  of  the  International  Conference  on  Lasers  *80,  New 
Orleans,  and  in  "Four-Wave  Mixing  Near  the  Band  Edge'',  which  is  submitted  for 
publication.  These  papers  are  reproduced  in  this  report  as  Appendices  1  and  2. 

We  examined  two-photon  resonances  for  the  purpose  of  four-wave  mixing 
and  found  that  the  third-order  susceptibility  is  typically  two  orders  of 
magnitude  less  than  what  is  practically  desired.  We  therefore  searched  for 
alternative  means,  particularly  for  phase  conjugation  at  the  iodine  laser 
wavelength,  1.3u. 

An  alternative  method  uses  the  interband  transitions  near  the  band 
edge.  We  found  that  the  reflectivity  for  conjugate  waves  is  resonantly  enhanced 
near  the  band  edge.  If  the  frequency  of  the  light  is  near  the  band  edge,  large, 
third-order  susceptibilities  and  reflectivities  are  obtained.  To  match  1.3u,  a 
particularaly  promising  material  turns  out  to  be  (Ga^  which  is  a 
semiconductor  compound  with  a  direct  band  gap  at  the  center  of  the  Bril  loin  zone. 
Th**  size  of  the  gap  can  be  varied  by  changing  the  alloy  composition.  It  is 
therefore  possible  to  match  the  band  gap  and  the  iodine  laser  frequency  as 
closely  as  desired.  The  reflectivity,  as  it  turns  out,  has  a  simple  Lorentz 
factor  near  the  band  edge,  which  leads  to  the  resonant  enhancement.  When 
X  -v  0.6,  the  (GajjIn^.jjlAs  band  gap  approximately  matches  1.3u.  The  third  order 
susceptibility  is  then  about  10”®  esu,  which  is  a  substantial  amount. 

The  calculations  we  have  performed  are  sufficiently  general  so  that 
they  can  be  adapted  to  other  semiconductors. 

We  also  considered  another  method  tr.  match  1.3u  with  a  single  photon 
resonance  of  a  nonlinear  medium.  Foreign  atoms  embedded  in  solids  or  liquids 
strongly  interact  with  their  surroundings  and  therefore,  the  relative  spacing  of 
their  atomic  levels  can  be  drastically  altered.  By  choosing  the  proper  solid  or 


liquid,  and  a  foreign  atom.  It  Is  possible  to  match  1.3y  with  an  atomic 

transition.  For  this  purpose  we  have  considered  alkali  Ions  embedded  In  solid 
dielectrics.  A  promising  system  like  this  has  an  additional  advantage.  Most  of 
the  absorption  at  1.3y  Is  due  to  the  Na^  themselves.  By  controlling  the  number 

of  Na^  In  the  sample.  It  1s  possible  to  have  a  system  which  can  easily  be 

saturated  at  l.Sy. 


III.  TRANSVERSE  MODE  COMPETITION  IN  GAIN  DRIVEN  UNSTABLE  RESONATORS 


The  analysis,  calculations,  results  and  discussion  are  presented  in 
Appendix  3,  which  is  a  preprint  of  a  paper  submitted  to  the  Journal  of  the 
Optical  Society  of  America.  Below  we  summarize  the  results  of  the  project. 

The  calculation  hinges  on  the  expansion  of  the  field  in  an  unstable 
resonator  with  a  saturable  gain  medium  in  terms  of  a  series  of  empty  resonator 
modes.  Substitution  into  the  Maxwell  equations  yields,  after  some  reductions,  a 
set  of  first  order  differential  equations  for  the  time  dependence  of  the 
expansion  coefficients.  The  diffraction  losses  are  included  in  a  natural  way  and 
the  results  are  physically  transparent. 

The  equations  are  coupled  together  by  overlap  integrals  involving  the 
empty  resonator  modes  and  the  nonlinear,  field  dependent  susceptibility.  The 
integrals  range  over  longitudinal  and  transverse  coordinates,  and  hence  are  quite 
complicated.  However,  if  the  empty  resonator  modes  are  assumed  to  be  known  in 
advance  and  can  be  stored  on  disc,  the  calculation  is  much  easier  than  the 
numerical  evaluation  of  the  nonlinear  Fresnel -Kirchoff  transform.  The  necessary 
accuracy  for  the  evaluation  was  determined,  and  it  was  found  that  the  grid  size 
requirements  for  the  longitudinal  integration  were  two  orders  of  magnitude  less 
stringent  than  for  the  transverse  integration. 

Calculations  were  carried  out  for  M  *  2  cavities  with  Ng^  *  6.5  and 
Neq  =  7.02.  These  two  values  of  equivalent  Fresnel  number  were  chosen  because 
they  are  the  locations  of  a  peak  and  a  crossing  point  of  the  eigenvalue  versus 
equivalent  Fresnel  number  plot.  The  lowest  loss  empty  resonator  modes  for  the 
Ngq  *  6.5  cavity  is  well  separated  from  the  others,  and  thus  we  expected  that 
this  mode  would  dominate  the  laser  operation.  In  the  calculations  this  indeed 
occurred,  with  about  95%  of  the  energy  concentrated  in  the  lowest  lots  modes. 
The  higher  order  modes  were  found  to  affect  strongly  the  output  coupling, 
frequency  pulling,  and  linewidth  of  the  laser,  but  not  the  distribution  of 
energy. 


By  contrast,  for  the  Ng^  =  7.02  cavity  the  two  lowest  loss  modes  have 
equal  losses,  and  the  posibility  of  unstable  behavior  was  considered.  Bistable 
operation  was  found  to  occur  in  a  small  region  of  pumping  strengths,  with  stable 
solutions  on  either  side  of  this  region.  In  the  region  of  small  gain  both  modes 
were  found  to  be  oscillatory,  and  above  the  bistable  regime  one  of  the  lowest 
loss  modes  was  suppressed.  In  the  region  of  bistability  itself,  self-sustaining 
pulsations  were  found  to  occur  for  most  initial  conditions  of  the  laser  system. 
The  iterative  solution  of  the  equations  of  motion  did  not  converge  to  a  stable 
point  in  these  cases,  but  to  a  well  defined  limit  cycle.  This  type  of  behavior 
is  well  known  in  laser  systems  and  is  a  result  of  nonlinearity.  This  is  the 
first  demonstration  of  nonlinear  oscillations  arising  from  competition  between 
transverse  modes. 

The  formulation  of  the  equations  of  motion  of  unstable  laser  resonators 
in  Lamb  Theory  form  has  been  shown  to  be  useful  both  qualitatively  and 
quantitatively.  Physically,  the  idea  of  expressing  the  laser  energy  distribution 
as  a  superposition  of  empty  resonator  modes  is  appealing.  Bistable  behavior  and 
self-sustaining  pulsations  result  naturally  from  this  model.  The  numerical 
results  are  expected  to  be  most  accurate  when  one  of  the  modes  manages  to 
suppress  all  the  others.  Quantities  such  as  frequency  pulling,  laser  linewidth, 
and  output  coupling  are  given,  as  well  as  the  stable  distribution  of  radiation  in 
the  cavity. 
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NOVEL  APPROACHES  TO  NONLINEAR  ADAPTIVE  OPTICS  AT  1.3u 


A.  £191  and  D.  Rogovin 

Science  Applications,  Inc. ,1200  Prospect  Street, La  Jolla,  CA  92038 


•Jbstract 


In  gases  degenerate  four-wave  mixing  at  1.3v  poses  unusual  difficulties  due  to  the 
lack  of  suitable  systems.  However,  nonlinear  adaptive  optics  at  these  wavelengths  appears 
to  be  quite  promising  in  solids. 


introduction 

Recent  advances  in  the  technology  of  the  chemically  pumped  atomic  iodine  laser  have 
generated  an  interest  in  the  resonantly  enhanced  phase  conjugation  via  degenerate  four- 
wave  mixing  at  1.3u.  Previous  studies  have  shown  that  the  reflectivity  of  the  conjugate  ... 
waves  is  sensitive  to  frequency  detuning  from  an  allowed  atomic  or  molecular  transition. 

For  instance,  the  reflectivity  for  NH,  is  essentially  zero  a  few  wavenumbers  from  an 
allowed  transition.  The  fact  that  one  needs  to  work  near  an  allowed  transition  limits 

drastically  the  number  and  type  of  materials  one  can  use  as  conjugators  at  1. 3u. Unfortun¬ 
ately,  there  do  not  appear  to  be  any  Icnown  gases  which  are  suitable  for  resonant 
degenerate  four-wave  mixing  at  1.3u,  at  least  not  in  thermal  equilibrium.  Althouch 
resonant  1 . 3u~transitions  do  occur  between  some  excited  states  in  atomic  gases,  their  short 
T.  times  as  well  as  the  additional  difficulties  due  to  maintaining  the  medium  in  a 
specified  nonequilibrium  state  make  their  use  impractical.  Two-photon  resonances  and  the 
use  of  neutral  iodine  are  also  impractical  as  discussed  below. However ,  there  are  solid- 
state  systems  which  are  quite  promising;  particularly  the  following  two: 

(1)  The  semiconductor  alloy  <Ga  In,_  )As  with  x  »  0.6, 

(2)  Atomic  or  molecular  impuritres''’”*embedded  in  solid  dielectrics. 

The  following  two  sections  are  devoted  to  the  discussion  of  these  two  systems.  In  the 
remainder  of  this  introduction  we  discuss  briefly  the  limitations  of  the  use  of  two-photon 
resonance  emd  neutral  atomic  iodine . 

To  examine  the  two-photon  resonant  enhancement  of  degenerate  four-wave  mixing  in  gases , 
we  consider  the  three-level  configurations  depicted  in  Fig.  1.  In  Confiouration  A,  levels 
0  and  2  are  not  coupled  by  radiative  transitions  and  uq  ^  ^  u  ,  but  2u  >  u.  .  If 

the  atom  or  the  molecule  is  initially  in  its  ground  state,  then  the  third-order  suscepti¬ 
bility  at  ut  is  approximately  given  by'*' 

X '  '  («)  »(3f5’)  Njuul  luul  )  (ui-uo+Iy)  *(ui-ui+iY  )  *  (2u-iuo-ii)  j +iY )  *.  (1) 

*  -1  (3) 

Typically/ Y  ^  2^10  rad. sec  at  1  torr  and  u  ^  u  I  Debye.  Let  (u  ^  1  eV  and 

2u  auo+ui  ,  then  at  Ohe  torr, 

I  X  (u)  I  ~  (6*10"^^esul  .  (u  +  0)  lu  -u  I  [  (u-  u  +4y^J"^^^  .  (2) 

0  1  0  1  0  1 


IM  ■) 


FIG. 1. Two-photon  resonance  in  three-level  systems. 


The  overall  ‘rsquency  factor  in  Eq. (2)  is  on  the  order  of . unitv  if  m  is  sianificantlv 
different  from  ^  .  To  achieve  useful  reflectivities,  I  should  b4  on, the  order  of” 

lij'^'‘  esu  or  ‘  larger.  To  increase  the  right  hand  side  of  (2)  to  lo"^'*^  esu  requires 
u /V  ui  <vuJo  <  in  which  case  single  photon  transitions  are  important  and  Eq.  (1)  does  not 
apply.  Besides,  in  that  case  one  has  found  a  system  whose  level  structure  matches  iodine. 

In  Configuration  B  ,  levels  0  and  2,  and  1  and  2  are  coupled  by  dipole  transitions; 

0  2md  1  are  not.  In  this  case  we  assume  that  2iu=  uo"  ui;  then,  at  one  torr, 

[6x10  ^^esul.(wo  “Oil  )^(  IJ0+  a>i)  (3) 

This  is  also  about  two  orders  of  magnitude  smaller  than  the  minimum  value  needed. 

The  impracticality  of  neutral  iodine  comes 
from  the  fact  that  the  lasing  levels  of  atomic 
iodine  are  coupled  by  magnetic  dipole 
transitions'  .  Since  magnetic  dipole 

transition  rates  are  smaller  than  electric  ,  j  | 

dipole  transition  rates  by  a  factor  of  (v/c)'^,  I  j 

the  corresponding  susceptibilities  are  related  by  1  / 

(v/c) Ig'tw)  •  (4)  \  r 

Since  (v/c)^10  ^  for  outer  electrons,  it  is  i \.  y'  ! 

clear  that  one  cannot  achieve  useful  reflectivities  j  i 

using  atomic  iodine.  i  [ 


use  of  (Ga jjln^_^) As 


When  w  matches  the  band  gap  of  a  direct  gap 
semiconductor,  the  four-wave  mixing  is  resonantly 
enhanced.  One  can  exploit  this  fact  for  phase 
conju'”tion  at  1.3u  using  (Ga  In,_  )As  alloys. 

They  are  direct  gap  materials^'and  the.fize  of  the 
band  gap  depends  on  the  composition  x.  '  To  a  good 
approximation  one  can  write 


E|.(x,T)=gg(T)*g^(T)x+g2(T)x'‘  .  (5)  |  j'' “ 

At  T=300K,  gg=0.36eV,  g.=0.77  eV  and  g,=0.31eV. 

It  is  therefore  possible  to  make  the  “^direct 

gap  as  close  to  1.3u  as  desired,  either  from  FIG. 2.  The  band  structure, 

above  or  from  below,  by  a  judicious  choice  of  ,,, 
the  composition,  and  thus  resonantly  enhance  x  • 

Further  fine  tuning  is  provided  by  temperature  since  o's  are  temperature  dependent, 

In  order  to  obtain  an  estimate  of  the  conjugate  wave  reflectivity,  we  calculate  v 
using  the  band  structure  shown  in  Fig.  2  for  (Ga  In._,.;As.  To  simplify  the  calculation  we 
assume  that  (a)  all  waves  are  polarized  in  the  5-dir^ction:  (b)  conduction  and  valence 
electrons  have  the  same  lifetime  y:  (c)  the  dipole  moment  is  given  by 


where  m„  is  the  bare  electron  mass,  c  is  the  hiah  frequency  dielectric  constant  and 
is  the  interband  momentum  operator:  (3)  intraband  transitions  are  negligible  and  only 
those  interband  transitions  which  are  near  the  band  gap  as  shown  in  Fig. 2  need  be  taken 
into  account:  and  finally  (e)  band  energies  are  parabolic: 


(E„/2)  -  (•n^k^/2m,)  ,  E  (k)  =- (E,,/2)  -  (-fi^k^ ,/2m  ) 

'J  C  V  {j  V 


for  v=l,2 
for  v=3 


is  given  by  (see  Ref.  2) 


(w)  =  (2-)‘^I  .'dk  !u^„!*(f(E  )-f(E  )1(4(L  -L  )('L  ^ 

V  cv  V  c  VC  cv  VC 


■)  +  (■f;.,*iv/2)  "■( 


where  »(«u-E  +E  ,+iY)“^,  f (E  )-(l+exp8(E  -E.)  1*^,  S-(k_T)"^,  and  E.  is  the  Fermi 

energy.  x  separated  into  a  t^mperScura-iadepmadent  and  a  tesi^rature-dependent 

part.  The  integrals  for  the  temperature- independent  part  can  be  evaluated  exactly. 
Designating  the  result  as  xA^'  >  on«  finds 

xP'  (u)-  XQn:j_/i_)”^[i+(2n  />(u)  +  (iY/a  )  +  (a  +iY)<>fw/2+iY/4)"^l 

-1  ~  1  <’) 
■^[l+2!5^(hu+iY)  ^+(iY/a^)  +  (i_-ix)  (Ku/2+iY/4) 

where  Xo*<27rh^)‘^  ^  I  ^ n^«  [  (Ej.ihu)  ^+Y^  ]  and  (E^+hui) . 

For  small  detuning  from  E^,  I E^-Huj /E^<<1,  the  +  terms  in  (9)  can  be  neglected  and 

Xo^'(w>  =  Xo(n_Ay^)"^(i.-iY).  (10) 

In  particular. 


!  xP'  (w) 


x'2xnip^^  -  /2xnrE_-hu+  -/(E.-Ku)  ^+y^] 


Thus,  lx  Q  I  rises  sharply  when  u  increases  above  the  band  gap  as  shown  in  Pig. 3,, 

The  temperature-dependent  part  of  the  susceptibility,  which  we  designate  as  xP  # 
cannot  be  integrated  explicitly,  although  an  approximate  form  for  |  E.-liu  |  <<E_  can  be 
obtained  readily.  At  thermal  equilibrium, E,  is  nearly  in  the  middle  of  the  gap  and  thus 
Ej  <<  Eg.  Using  this  fact,  one  finds  ^ 

xP’(u)ti  -d+i)  lu^|S2m^m^)^/2(m^+m^)‘^/2lexpt-S(Eg/2-Ej)]  (12, 

"  y®  (  (y^+Eg-hu)  ^♦Y^l  ^{exp  (-Sm^y^  (m^+m^)  ”^] +exp  (m^+m^)  ) 

IS (1+i) (k„T) (e"® (Eq/Z-E^ ) _g-3 (Eg/2+E ,) 

lOp - 3  =1 - 2 - 

:  <v/2v»-1  ,v(3),  :  4P/VYt(Eg-Kw)WP 


In  obtainina  (12)  ,  we  neolected  (Eg-Wu)  or 
Y  which  appear  toget.her  with  y‘  in ’the 
numerator,  and  also  used  the  fact  that  v<< 
Eg  or  Hu.  It  is  clear  from  (12)  that 

lxP’l/!xP’-(V/^p/V'®V2>  ,  ,13, 


is  extremely  small  compared 


-6.0  -4.0  -2.0 


4.0  6.0 


FIG. 3.  The  third-order  susceptibility 
as  a  function  of  the  frequency 
detuning  ■;=(£  -hw)/-.-. 


to  Xn  due  to  the  e.Yppnential  factor.  For 
all  intents,  x'^^  =  Xp^.  ,3, 

We  can  now  estimate  x  .  Det  us  consider 
a  particular  comoosition  with  x-0.6  at 
T=300K,  for  which  Eg=0.95  ev  and  thus  matches 
1.3u.  The  other  parineters  arec  =10,  m.  =  0.1mg, 
nu^O.Sm.  andiu  i  =10  Debye.  Tne  finite 
lifetime  for  electrons  and  holes  arises 
mainly  from  phonon  scatterings  and  is  on  the 
order  of  /»*  10  -10  sec.  Using 

these  values  one  finds  that  at  1.3u, 
lx'  I<wl0~  -  4*10  esu.  These  are  extremely 
promising  values. 

Since  both  object  and  conjugate  waves  are 
damped  in  a  lossy  medium,  the  reflectivity 
also  depends  on  absorption  even  when  Dump 
waves  are  assumed  to  ba.undamped.  The  ref¬ 
lectivity  is  aiven  by 


0  ■  licsinai^ecosa  ♦  oslncif^ 


k|  •  16,r2<^*2n-2^|^(3)|^  (X4) 

5  •  -  o}, 

where  t  is  the  length  of  the  nonlinear  mixing  region,  is  the  index  of  the  refraction, 

I.  and  I.  are  the  intensities  of  the  two  counterpropagating  pump  waves,  a  is  the  absorption 
coefficient.  At  room  temperatures  a  can  be  written  approximately  as  '  ' 


5  “0"“l 


-Oj^  (Eg-Ku) 


for  Ku<Eq 
for  )<ii)>E_ 


where  a„al.3i<10^  cm~^ 


a,  =10^  cm”^eV  ^  and  a«»6>«10^  cm~^  ev”^^^  .  It  is  clear  that 


absorption  will  be  high  and  one  willgusually'^be  in  the  absontive  regime,  i.e..C  will 
usually.be  imaginary.  Let  |x'  M^IO  esu  and  MH/cm^ ,  then  |ic|’'260  cm  .a  is  a 
10^  cm  .  !Vhen  Sis  imaginary  n  is  independent  of  c  and  is  given  by 


Ti=|«:|^/(2a)^  =  2x10"^-  2% 

This  corresponds  to  a  reasonably  efficient  phase  conjuaation. 


impurities  embedded  in  solid  dielectrics 

It  is  clear  from  the  above  discussion  that  if  thtir  linewidths  were  smaller,  solids  would 
be  extremely  efficient  for  fpurg^ve  mixing.  One  approach  for  achieving  narrow  linewidths 
in  solids  is  to  use  discrete^X/hidn  couple  weakly  to  phonons  compared  to  continuum  states 
of  the  energy  bands.  For  exaunole,  when  an  atom  or  an  ion  is  embedded_in  a  dielectric, 
the  lifetimes  of  the  low  lyino  states  are  usually  on  the  order  of  10~  sec.  One  can  also 
use  the  fact  the  atomic  energy  levels  are  reduced  approximately  by  a  factor  cT  ,  where 
e.  is  the  static  dielectric  constant,  to  tune  some  atomic  transition  frequencies  to  1.3u 
by  choosing  e«  appropriately. 

For  the  sedte  of  concreteness  we  consider  alkali  type  ions  embedded  in  semiconductors 
and  insulators  in  the  followina.  However,  the  basic  idea  can  be  extended  to  neutral  atoms 
and  other  types  of  ions  in  hichlv  degenerate  solids  or  glasses.  Also,  the  system  whose 
discrete  states  are  to  be  used  need  not  be  an  ion  or  atom;  it  could  equally  well  bq  a 
defect  in  the  solid. 

The  discrete  energy  levels  of  an  alkali  ion  can  be  approximated  by' 


E„  -  E_  -  R(Z-Z'+l)^er^(n-i)"^  (17) 

n  G  w 

where  R  is  the  Rydberg  enercy,  Z  is  the  nuclear  charge,  Z'  is  the  number  of  electrons, 
n  is  an  integer  and  a  is  the  quantum  defect.  The  energy  levels  are  measured  from  the 
edge  of  the  valence  bands.  For  resonant  tuning  we  require  1. 3u“K«*0. 9SeV»E  -E  ,  for 
some  n' >n.  Thus, 


(Z-Z’+l)^[(n-a)'^ 


(n'-a)"^]  . 


One  way  to  proceed  is  to  choose  a  specific  atom  or  ion,  as  well  as  a  particular  transition 
and  then  solve  for 

Eg  =  3.8(Z-Z’-t-l)  ln'-ar^tn-ar^(n’-n)^''^(n+n'-2a)^^^  .  (191 


Below  we  list  the  values  of  £-  given  by  Eo.  (19)  for  a  number  of  ions. 
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s  al  1° 

Na'^:  Z-Z'+l-2,  i-1.4  1  2  12 

1- 19 

2  3  12 

2- 13 

Al*'*'*;  2-Z'+l-4,  i-0.9  1  2  68 

1  -  152 


P®"^:  Z-Z'+l-6,  a-0.7  1  2  99 

1-76 


As  static  dielectric  constants  qo,  all  of  these  are  reasonable  values.  In  small  gap  mater 
ials  (E_'^  0.2eV)  c.  is  usually  large.  For  insulator  and  hiah  gap  materials  e.  is  on 
the  ordef  of  10.  ^  “ 

At  first  glance,  it  may  seem  that  there  would  be  strona  absorption  in  small  gap  materi¬ 
als  like  PbSe  or  PbTe  alloys,  since  |4u>C.  and  thus  they  would  have  negligible  reflec 
tivity.  However.  leV  is  so  large  compared  to  the  gaps  of  these  materials  that  the  tran¬ 
sition  would  fall  onto  very  high  conduction  bands.  Bv  a  slight  adjustment  of  the  composi¬ 
tion,  the  relative  position  of  these  bands  can  be  altered  and  absorption  decreased.  Cq 
can  also  be  adjusted  by  varyinc  the  alloy  composition. 

The  numbers  for  Na*  are  quite  close  to  pure  Si  and  CdS.  Si  has  e«»11.8.  CdS 
is  an  insulator  and  has  E_>2.4eV.  CdS  has  an  added  advantage  in  that  absorption  at  1.3u 
would  only  be  due  to  the  iiipurity  sites  themselves,  and  hence  these  transitions  can  be 
saturated  for  sufficiently  low  impurity  densities. 

We  acknowledge  the  support  of  AFOSR  for  this  work. 
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ABSTRACT 

We  discuss  degenerate  four-wave  mixing  and  phase  conjugation 
near  the  band  edge  of  direct  gap  semiconductors  and  show  that 
the  reflectivity  for  conjugate  waves  is  resonantly  enhanced 
in  the  vicinity  of  the  band  edge.  We  also  compare  the  quantum 
mechanical  results  with  the  classical  Orude  model. 


I.  INTRODUCTION 


There  is  presently  a  considerable  body  of  work  on  phase  conjugation 


via  degenerate  four-wave  mixing  in  semiconductors. 


(1-6) 


It  has  been  moti¬ 


vated  partly  by  the  availability  of  a  variety  of  frequency  ranges,  and  partly 

by  the  relatively  high  efficiencies  observed  in  these  materials.  This  paper, 

too,  is  concerned  with  degenerate  four-wave  mixing.  We  consider  a  direct 

gap  semiconductor  and  calculate  the  contribution  of  interband  transitions 
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near  the  band  edge  to  the  third  order  susceptibility  x  v  (XJ  ®  UJ  +  (J  -  U))  . 

We  also  consider  the  absorption  coefficient  a  in  the  semiconductor  in  order 
to  obtain  a  realistic  expression  for  the  conjugate  wave  reflectivity  n. 


There  are  two  reasons  for  the  interest  in  a  band-edge  mechanism 
for  four-wave  mixing.  First,  is  resonantly  enhanced  as  Tioi  increases 
i;o  Eg,  the  band  gap.  As  a  result,  n  is  also  enhanced  despite  increased 
absorption  (which  adversely  affects  phase  conjugation  ’).  Second,  a  band- 
edge  mechanism  offers  a  greater  potential  for  frequency  tunability,  since  in 
many  semiconductor  compounds  such  as  Ga,  (InAs)  ,  Cd,  (HgTe)  ,  Pb,  (SnSe)  , 

X  A*X  X  A*X  X 


etc.,  Eg  can  be  continuously  varied  by  changing  composition 


(9) 


A  number  of  authors  have  discussed  band-edge  mechanisms  for  four- 
wave  mixing. ^ Among  these,  R.  K.  Jain,  et  al.,^^’^^  have  considered 
a  mechanism  which  is  similar  to  what  we  have  in  mind.  They  imagine  creation 
of  an  electron-hole  plasma  in  the  semiconductor  and  use  a  Drude  model  to 
calculate  In  this  paper,  however,  we  perform  a  density  matrix  cal¬ 

culation  to  provide  a  more  precise  pictu’^e  of  the  frequency  dependence  of 
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» 


♦ 


f 


t 


(3) 

X'  and  n.  Also,  our  results  apply  to  frequencies  that  are  lower  than  E-. 

b 

A  detailed  comparison  of  our  results  with  those  of  Ref.  1  is  given  in  Section  V. 

Finally,  we  also  note  that  M.  A.  Khan,  et  al.,^^^  have  studied  a  case  in  which 
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X'  is  resonant  for  ■huj  Eg  because  of  a  nonparabolic  band  structure.  This 
nonlinear  process  is  quite  different  than  the  one  taken  up  here.^^^^  Here, 
we  are  concerned  with  an  ordinary  and  direct  electron-field  coupling  via 
interband  transitions. 

In  Section  II,  we  describe  the  band  structure  used  in  the  calculation. 
It  is  a  simplified  version  of  the  band  structure  of  III-V  compounds,  chosen  with 
the  phase  conjugation  of  iodine  laser  beams  in  mind.  In  Section  III  we  give 
the  general  results  for  a  two-component  plasma  in  the  semiconductor.  In  Sec¬ 
tion  IV,  the  general  results  are  specialized  to  small  perturbations  from 
thermal  equilibrium.  Finally,  in  Section  VI,  we  discuss  the  reflectivity  of 
a  spacific  direct  gap  semiconductor,  Ga^_^( InAs)^, 

The  main  result  of  our  work  is  that  the  reflectivity  of  conjugate 
waves  has  a  resonance  in  the  vicinity  of  the  band  edge.  The  precise  position 
of  the  resonance  depends  parametrically  on  the  sample  thickness.  The  resonance 
shifts  to  lower  frequencies  for  thicker  samples,  as  one  might  expect  from  the 
depletion  of  the  pump  waves.  If  the  depletion  of  the  pump  waves  are  negligible, 
as  when  the  semiconductor  is  saturated,  then  the  reflectivity  is  nearly  constant 
above  the  band  edge  and  decreases  rapidly  below  it,  independently  of  the  sample 
thickness. 


» 
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II.  BAND  MODEL 


To  be  definite,  we  assume  a  band  structure  as  shown  in  Fig.  1, 
which  is  typical  of  III-V  compounds.  It  has  one  light  hole  band,  two 
heavy  hole  bands  which  are  degenerate,  and  one  conduction  band.  To 
simplify  the  calculation  we  also  assume  that  the  bands  are  parabolic; 

E^(k)  =  (Eq/2)  +  (-n^k^/am^M 

Ey(k)  =  -  (Eg/2)  -  (liV/2m^*),  (v=l,2.3)  . 

where  the  m*'s  are  the  appropriate  effective  masses  near  the  center  of  the 
Brillouin  Zone.  For  the  calculation  one  also  needs  the  dipole  matrix  element 
for  the  interband  transitions.  It  is  related  to  the  more  familiar  interband 
momentum  matrix  element  by  the  following  relations: 


Here  is  the  high  frequency  dielectric  constant,  m^  is  the  bare  elect'^on 

-1/2 

mass,  and  is  the  one-electron  crystal  Hamiltonian.  The  factor  c 
0  * 

in  (2a)  comes  from  high  frequency  screening  of  electrons  and  renormalizes 
the  bare  electronic  charge.  Since  our  concern  is  with  the  band  edge,  the 


right  hand  side  of  (2c)  need  be  evaluated  only  at  k=0.  can  be  taken 

as  constant. 

We  are  interested  in  those  situations  in  which  the  semiconductor 
is  either  at  room  temperature  or  at  a  temperature  which  is  not  too  low 
relative  to  room  temperature.  At  these  temperatures  the  absorption  edge 
is  smoothly  broadened  to  frequencies  which  are  considerably  below  the 
band  edge,  due  to  interactions  among  electrons,  holes  and  phonons.  To 
take  into  account  this  broadening,  we  assign  average  linewidths  and 
to  electron  and  hole  states,  respectively.  This  procedure  should  be 
adequate  to  describe  the  actual  situation  for  intermediate  and  room  tem¬ 
peratures.  However,  for  lower  temperatures,  there  are  well  defined  discrete 
excitonic  peaks  and  one  must  take  the  discrete  nature  of  these  peaks  into 
account.  This  requires  considerable  modification  of  the  calculation  and 
is  reserved  for  another  publication.  Here  we  confine  ourselves  to  the 
regime  in  which  discrete  excitonic  peaks  are  not  observable. 
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III,  THIRD  ORDER  SUSCEPTIBILITY,  ABSORPTION  AND  REFLECTIVITY 


I 


I 


» 


» 


I 


) 


The  geometry  of  degenerate  four-wave  mixing  is  shown  in  Fig.  2. 

We  assume  that  all  four  waves  are  linearly  polarized  in  the  z-direction 

and  actually  calculate  X^^zz  •  Since  absorption  is  important 

for  phase  conjugation,  we  also  calculate  the  frequency-dependent  absorption 

coefficient  a(a))  under  the  same  assumptions.  A  density  matrix  calculation 
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gives  the  following  expressions  for  and  ct  {see  Ref.  8  and  Footnote  11): 


.(3)  . 


Z  I  dk(2.)-3j4|^[/(E^)  -  fiE^)] 

X  [(Yc+Yv)^(YcYv)'^(Lvc-‘-cv^^i'-vcl^  ^  |L^J^) 

+  (41iu)+iY(.+iYy)(2h6a+iY^)‘^(2riw+iYy)"^(L^y-L^^)J  , 


(3a) 


a  =  W(cTr  E  )‘-^  Y.  Ycvl^cv 

^  v=l  BZ 


V  V  I.,2  |2  /dirr>(E  )-f(E  )] 

5Z  L  V  c  j 

X  [l^cvl'  ^  l^c!'] 


where  n  =  /I  is  the  index  of  refraction,  and 


^nn'  =  ^  ^  ’'^nn']’^  ' 

Ynn^  =  (V^n'^/2  , 

;(E^)  =  [l  +  exp8(E^(i?)-E^^)]-^  .  B  =  (kgT)-^  , 


Efn 


E,  for  n=c 
fc 

-E^^  for  n=v=l,2,3  . 


(3b) 


(4a) 

(4b) 

(4c) 


(4d) 
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In  (4c, d),  are  Fermi  energies  for  electrons  and  holes,  assumed  to  be 
distinct  in  order  to  include  those  situations  in  which  absorption  at 
significantly  affects  band  populations. 

One  can  let  the  upper  limit  of  k  go  to  infinity  to  simplify  the 

integrals.  We  do  not  expect  much  error  resulting  from  this  procedure  since 

as  k-^“,  the  integrands  decrease  rapidly  for  parabolic  bands.  However,  the 

situation  is  quite  different  for  nonparabolic  bands.  For  instance,  the 
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integrals  for  Kane's  two-band  model  ’  diverge  if  one  lets  the  upper  limit 
of  k  go  to  infinity;  therefore,  they  must  be  evaluated  for  a  finite 
Brillouin  Zone. 

X  and  a  consist  of  two  types  of  terms, 

=  L  (K°  +  kJ),  (5a) 

3 

a  =  E  (A°  +  AM,  (5b) 

v=l 

where  and  A°  are  temperature  independent,  and  and  A^  are  temperature 
dependent.  This  separation  is  effected  by  writing  the  population  factor  as  ’  - 

/(E^)  -  /-(E^)  =  1  -  Cl-f(E^)  +  /(E(.)]  (6) 

=  1  -  [[l  +  exp3  (-E^-E^^)]‘^+[l+expS(E^-E.^)]'^  j  . 

The  integrals  for  K°  and  A°  can  be  evaluated  by  contour  integration  (see 
the  appendix)  and  one  finds 
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(7a) 


+  [(V>cv^(^+^^’^]tL(Mjv)'^-(^^‘^'^2iY^.)‘^-(41la)+2iYy)’^] 

-  [2^cv(Vv^‘"’^^'^]l 


and 


where 


and 


aJ  =  4na32^J/2|^Zj2(^2^^^j-lJ-^^^2(-nu3-Eg)] 


1/2 


+  (Egirhoj) 


(7b) 


(7c) 

(7d) 


(7e) 


The  temperature-dependent  terms  cannot  be  evaluated  explicitly 
for  arbitrary  temperature.  We  write  them  in  terms  of  Fermi -Dirac  integrals 


of  the  dimensionless  form 

fJJ  (j;a)  =  j  dy  y"(y^+a)""’[^l+Cj  exp  (bjy^)J  .  (8) 

Let  us  define 

^  =  (EQ-*-iY^y)Ahw  ,  (9a) 

c^  =  expS  (Eg/2-E^^)  ,  (9b) 

C2  =  exp  3  (Eg/2-E^^)  ,  (9c) 

bi  =  (m^m*)  STiai  ,  (9d) 

^2  =  ("’v/'''v)  37iuj  .  (9e) 
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Ky  and  are  then  given  by 


2 

X  Zjlm[F2(j;^-l)  +  {j;i;*+l)] 

+  i(Y^/M[F2(j;^*-l)  -  F^(j;c+1)] 

+  (Yj;^/'hu))Re[i(c-l)F°{j;c-l)  +  i  (c*+l)F  j(  j  ;c*+l)] 
+  (Yjy2hu))[(2Tic.+  iY^)“^  +  (2tico+iYj'^][F^{j;!;+l) 


;(j;;*-i)]] 


(lOa) 


=  4(jjn  I  1  ^ 

V  '  cv ' 


y 

(2m^cu/?l)^/^(7;cEg)'^  +  F^(j ;;+!)] 


Finally,  the  reflectivity  for  conjugate  waves  in  degenerate 
(7  31 

TTiTvifin  TC  ni\/nn  K\/'  • 


four-wave  mixing  is  given  by' 


where 


2  -2 
=  |<sinEx.l  lEcosEi  +  a  sinC?.|"  , 


|<|  =  16.2.c-V^(t  T  )1/2,y(3) 


(11a) 


(lib) 


^  -  i\  |2  2a/2 

5  =  (  1<|  '  3t  )  ' 


(lie) 


In  Eq.  (11a),  %  is  the  length  of  the  nonlinear  mixing  region,  and  and 
I2  are  the  intensities  of  the  two  counterpropagating  pump  waves.  We  note 
that  (11a)  is  derived  under  the  assumption  that  the  conjugate  and  signal 
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waves  are  depleted  but  the  pump  waves  are  not.  When  the  depletion  of  the 
pump  waves  is  also  taken  into  account,  as  seen  from  Fig.  2,  the  constant 
pump  amplitudes  and  A2  are  replaced  by 

A^(x)  =  A^(£)  exp  [-(£-x)a/2]  , 

A2(x)  =  A2(o)  exp  (-ax/2)  , 
and  thus  the  product 

A^(x)A2(x)  =  exp  [-al/2] 

is  independent  of  position,  [k]  becomes 

!<|  =  16TT2^c'^£;^(Ijl2)^''^iX^^^exp  (-a£/2)  . 

4  -1 

Above  the  band  edge  a  is  usually  on  the  order  of  10  cm  .  All  four 
waves  can  be  severely  depleted  and  one  should  use  (llg).  The  use  of 
(lib)  should  be  restricted  to  frequencies  that  are  lower  than  the  band 
edge.  There  is  still  significant  absorption  at  these  frequencies, 
however,  due  to  the  broadening  of  the  absorption  edge. 


(lid) 

(lie) 

(llf) 

(llg) 
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IV.  MODERATE  EXCITATIONS 


In  order  to  obtain  efficient  phase  conjugation,  one  would 
naturally  seek  a  situation  in  which  absorption  at  is  as  low  as  possible. 
We  now  take  this  to  be  the  case  and  assume  that  electron  and  hole  popula¬ 
tions  are  not  perturbed  significantly  from  their  equilibrium  forms.  We 
note  in  passing  that  this  assumption  contrasts  with  the  premise  of  Ref.  1 
in  which  Jain  and  Klein  imagine  creation  of  a  dense  electron-hole  plasma 
in  the  semiconductor.  In  their  case,  it  would  be  more  appropriate  to  assign 
a  separate  Fermi  energy  to  each  component  of  the  plasma,  since  each  ther¬ 
mal  izes  separately  before  the  interband  transitions  can  take  over  and  bring 
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the  two  components  into  a  common  thermal  equilibrium.'  ' 


At  the  thermal  equilibrium,  there  is  one  common  Fermi  energy: 


'f  =  ^fc  =  -^fh 


is  nearly  in  the  middle  of  the  band  gap,  given  by 


=  (k3T/2)tn[l*E  . 


One  can  therefore  set 

c^  =  C2  =  exp  (SEg/2)  »1 

for  band  gaps  on  the  order  of  0.5-1  eV  and  for  room  temperatures.  Also, 
the  integrand  in  (8)  contributes  mostly  in  the  vicinity  of  y=0  and 

F"(j;a)=  (2Cj)'^a-%j’(V'^  . 


(12b) 


.  ^  =  y  and  focus  our  attention  on  the  vicin 
which  we  mean 

the  vicinity  of  the  band  ed9e  (7a)  yields 

3/2  5  -1/2 


(15a) 


where 


5  =  (E--1iw)/x 


(15b) 


and 


X  =  (2ttTi  Y 
0 


V'^)-‘  i:  "’ri“cv'‘‘  • 


(15c) 


From  (7b),  one  has 

£  A”  .  • 


(163) 


where 


'^cv' 


=  4nu)‘^Y  ^ 


(16b) 


Similarly,  the  use 


•  nf  the  band  edge  yields 

of  (14)  in  the  vicimty  or  rne 


Ty}  =  Xa.(5+l)(5+i)(l-^^-^^  ^ 


(17a) 


and 


£a)  = 


(17b) 
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> 


I 


where 


«  2(i.l)(2iT6)'^''W^€xp(-6Eg/2)i;iu^y!^[(in*)^''2+(m*)^/2]  (17c) 


We  can  now  compare  temperature-dependent  and  temperature-independent 


terms.  From  (15c),  (16b),  (17c)  and  (17d),  one  finds 


IXt/Xol  (2^)(k3T/7rY)^'^^  exp  {-BEg/2) 


|at/aol  (2TT)(1ia)/2T7Y)^^^  exp  (-6Eg/2) 


{18a) 


(18b) 


Due  to  the  exponential  terms,  the  temperature-dependent  terms  are  clearly 

negligible  either  for  room  temperatures  under  thermal  equilibrium  conditions, 

or  for  small  perturbations  from  thermal  equilibrium.  For  lower  temperatures 

(31 

they  decrease  further,  exponentially.  For  all  intents  and  purposes,  x  and 
a  are  given  by  the  temperature-independent  terms  under  moderate  excitations 
as  long  as  the  band  gap  is  not  too  small. 


Neglecting  the  temperature-dependent  terms. 


> 


t 


t 


» 


f 


I 


» 


a  .  a I/U?  - 

O'*  -* 

=  +  <5]"^^^  . 

Rex^^^  =  {Xq/%)(1+<5^)"^'^^  ct 
0  3r  4  ^  Al-l 

■  V  L%  ^  “  J  ’ 

and 

=  (Xo/ao)(l+6^)^/^  ct^ 

=  (2xQ/aQ)a^[aJ  +  . 

Finally  one  needs  for  the  reflectivity: 

Ix^^^l  =  ^  Xo(l+5^)'^'''^(/l^+  6)"^  • 


(19a) 

(19b) 

(19c) 

(19d) 

(19e) 

(19f) 


(19g) 


Fig,  3  shows  |x^^^  1/(‘^Xq)  3S  a  function  of  6  =  (Eg-'Fiu))Y'^,  the 
frequency  detuning  per  linewidth.  Below  the  band  gap,  6»1  and 

ix^^h  =  2Xq(26)-^/2  ,  (20a) 

which  should  be  compared  with 

a  =  ao(26)"^  (20b) 


I 
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in  the  same  region.  Thus,  decreases  faster  compared  to  a  for 

(S>  1.  Above  the  band  gap,  <$<  o,  |6|»1,  and 

=  2Xo(26)^/2  ^ 

which  has  the  same  6-dependence  as  a. 


I 


(20c) 


I 


V.  COMPARISON  WITH  DRUDE  MODEL 


We  can  compare  our  expression  for  with  that  of  Jain  and 
Klein^^^.  Their  expression  for  obtained  from  a  Drude  model,  is 

X^Q^=  nce^Tn*a*(8TTm*^-huj^)"^  (21) 


Here  t  is  the  lifetime  for  electrons  and  holes,  a*  is  the  total  absorption 
coefficient,  {n*a*)  designates  the  fraction  of  the  absorption  which  corre¬ 
sponds  to  an  actual  electron-hole  pair  creation  (i.e.,  ri*a*  corresponds  to 

our  a),  and  m|^  is  the  reduced  mass  of  an  electron  and  a  hole;  m*^  = 

-1  ( 3l 

m|  mjij  {m*  +  m*)  .  Xp  is  real  and  describes  an  adiabatic  response  of 

the  system.  It  corresponds  to  our  Rex^^^.  Our  expression  for  has 

(31 

an  additional  part  which  is  imaginary.  Im  x  describes  those  processes  in 
which  at  least  one  real  transition  occurs.  Such  processes  can  either  put 
energy  into  the  system  or  extract  energy  from  the  system.  Note  that  above 
the  band  edge,  6<0  and  l6l>>l.  One  therefore  has  and  (19f) 

simplifies  to 


iTix^^’  =  (2Xo/aQ)cx 


(22) 


( 3l 

Thus,  above  the  band  edge  Imx'  '  is  proportional  to  the  absorption  coef¬ 
ficient,  having  the  same  functional  dependence  on  the  frequency  detuning 
(Eg-Tioi)  as  a.  This  is  consistent  with  having  at  least  one  real  transition 
occur. 
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To  compare  and  Rex^^^,  we  make  the  following  approximations. 
The  ratio  of  the  sums  over  the  dipole  moments  can  be  estimated  as 


Using  (2c) ,  one  has 


cl' 


€V|P^j|2(m2nV)'*  • 


The  momentum  matrix  element  is  related  to  the  electron  and  hole  masses 
near  the  band  edge  (see  Ref.  12); 

With  these  approximations  and  also  using  y  -  (fi/t),  one  finds 
Rex^^^  «  (■Fiu)/2)[(-nu-Eg)^  +  ’ 

Thus,  as  far  as  the  frequency  dependence  is  concerned,  the  two  expressions 
differ  only  by  the  square  root  of  a  Lorentz  factor.  The  other  factor  n”^ 
arises  from  two  sources.  One  is  that  we  are  taking  into  account  the  short 
wavelength  screening  and  treating  electrons  as  particles  each  wi ch  an 
effective  charge  =  e/n,  as  in  (2a).  This  contributes  a  factor  n 

in  (25).  The  remaining  n  factor  comes  from  the  difference  between  the 
way  we  count  the  number  of  photons  and  the  way  it  is  done  in  Ref.  1.  We 


use  the  standard  expansion  for  the  electric  field  operator  such  that 


(see  Ref.  14) , 

?(x,t)  e^|bj^exp^i(tt^*x-u)j,t)J>b^exp|ji(]<jrx-uj^t)J| ,  (26a) 

with 


where  X  designates  the  modes,  e^  is  the  unit  polarization  vector,  V  is 
the  volume,  and  and  are  related  by 

01^  *  (c/n)k^  .  (26c) 


For  a  single  coherent  mode  at  w, 

(4TT)'^/dx<E2^=  tlw<b^b>  »  tiwN  , 

where  N  is  the  average  number  of  photons  in  this  mode.  Thus,  the  photon 
current  density  is 

i 

t 

;  For  complex  field  amplitudes, 

•*  *  ' 

1  \ 


This  expression  for  differs  from  the  corresponding  one  in  Ref.  1  in 

-2  *4  (3) 

that  (29)  has  an  extra  n  factor.  Finally,  if  we  write  out  n  Xq  » 


-4 

we  see  that  n  factor  helps  to  renormalize  both  the  charge  and  the 
speed  of  light. 

It  is  interesting  that  the  presence  of  n"^  somewhat  compensates 

for  the  Lorentz  factor  in  (25)  and  makes  Rex^^^  and  comparable,  since 
2 

usually  n  '\-  10-20  for  semiconductors. 
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VI.  REFLECTIVITY  FOR 

We  now  consider  a  specific  semiconductor;  Gaj_^(InAs)^  with 
X  =  0.4  for  which  the  band  gap  is  Eg  =  0.95  eV.  This  coice  is  motivated 
by  the  fact  that  =  ETrticEg^  =  1.3  ym,  which  coincides  with  the  wave¬ 
length  of  high  energy  iodine  lasers.  There  has  been  a  practical  interest 
in  obtaining  phase  conjugation  at  the  iodine  wavelength.  To  the  authors' 
knowledge  there  is  no  gaseous  or  atomic  system  which  matches  this  wave¬ 
length  and  which  might  be  used  for  resonant  enhancement  of  the  conjugation 
process.  Neutral  iodine  itself  cannot  be  used,  since  the  lasing  transition 
is  a  magnetic  dipole  transition.  A  semiconductor  such  as  Ga^_j^( InAs)^  is 
therefore  a  natural  choice  for  phase  conjugation  at  X^  =  1.3  ym. 

The  remaining  parameters  for  our  sample  are  =  10,  m*  =  m^  = 

0.1  and  m*  =  m^  =  0.6  m^  (see  Ref.  15).  A  typical  lifetime  is  t  =  (Ti/y) 
s  10*^^  sec'^.  Thus, 

m*.  =  m*m*(m*+rnt) =  0.09  m^ 
eh  c  1  c  1  0 

and 

=  etl(2m*^£^  Eg)"^  =  10  Oebye, 

where  we  used  (23b)  and  (24).  Setting  Eq  =  "huj  in  (15c)  and  (16b)  yields 
-A  3  -1 

¥  *  3  X  10*  esu  and  :i  =  1.2  x  10  cm  .  The  value  of  a.,  which  corre- 
^0  0 

spends  to  the  value  of  ct  at  the  exact  band  edge  i  =  (Eg-?iu;)/Y=0,  is  in 
reasonable  agreement  with  the  experimentally  observed  value  for  GaAs  (see 


(31a) 


(31b) 
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page  61  of  Ref.  15).  Further,  let  I j=l2=lMW/cm^,  then  from  (lib),  |k1» 

1.1  X  10^  cm  which  is  comparable  to  but  less  than  a.  Two  comments  are 
in  order.  The  first  is  that  jx  and  hence  !<!  are  quite  sensitive  to 
the  linewidth.  For  example,  for  t  =  10'^^  sec'^,  one  finds  a  4  x  10^  cm"^, 
Xjj  «  10  ^  esu  and  1<1  »  36  cm”^.  The  second  is  that  the  parameter  ^  of  (11c), 
which  actually  determines  the  reflectivity,  is  usually  an  imaginary  quantity. 
For  imaginary  5,  (11a)  simplifies  to 


n  =  lKl^(2a)'^  . 


Thus,  if  the  pump  waves  are  not  depleted  and  x  =  10'^^  sec”^,  one  has  from 
(32)  and  (11b)  that  n  =  0.22  =  22%,  which  is  a  reasonably  efficient  phase 
conjugation. 


The  frequency-detuning  dependence  of  n  as  given  by  (32)  can  be 
written  out  as 


n 


exp 


/~7 


where 


Hq  =  27t  Eg  (t  u)  c  Y  n  ) 


2  ,2_2„4_6x-l,  v'_3/2|  2  , 4,2,  ^  „3/2  ■ .  z 


m 


'^cv ' 


cv 


12 


and 


We  used  (llg)  to  obtain  (33a).  The  parameter  may  be  called  the  optical 
thickness  of  the  sample.  If  one  wants  to  ignore  the  pump  depletion,  one  can 


(33a) 


(33b) 


(33c) 
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let  Pq  ->•  0  in  (33a),  which  would  then  revert  to  the  form  that  is  obtained 
from  (11b).  Figure  4  shows  (h/nQ)  vs  6  for  various  p^.  Above  the  band 
edge, 

(n,  Hp)  =  2Q-(4(52)“13exp  (-p^/2|6|). 

Thus,  when  the  pump  waves  are  undepleted,  n  has  an  upper  limit  which  is 
ZHq  for  frequency  detunings  that  are  large  relative  to  the  linewidth.  For 
the  specific  semiconductor  above  and  for  =  I2  =  1  MW/cm  and  t  =  IC''^’  sec, 
2no  =  .44  =  44%.  Similarly,  below  the  band  edge, 

(h/Hq)  =  (26^)"^  exp  {-oJv'zE), 

which  decreases  rapidly  with  increasing  5. 

It  is  clear  that  the  optimum  place  for  phase  conjugation  is  in 
the  vicinity  of  the  band  edge.  When  the  pump  depletion  is  taken  into 
account,  there  is  more  structure  in  the  reflectivity  as  a  function  of  the 
frequency  detuning.  Also,  the  reflectivity  decreases  as  tiu)  -  Eg  increases 
above  the  band  edge.  The  reflectivity  generally  peaks  in  the  vicinity  of 
the  band  edge.  The  position  of  the  peak  is  given  by  a  fifth  order  equation 
which  is  obtained  by  setting  the  derivative  of  (’^/’^q)  with  respect  to  5 
equal  to  zero; 

y^  +  y  -  (4/0^)  =  0 

where 


(34a) 


(34b) 


(35a) 


(35b) 
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There  is  at  least  one  root  of  Eq.  {SaS),  say  y^,  which  is  real  and 
positive.  The  left-hand  side  can  be  expanded  as 

y^  +  y  -  (4/Pq)  =  (y-yQ)g(y)  ,  (36a 

g{y)  =  y'^  +  y^y^  +  +  y^y  +  y^^^  +  l  -  (36b 

The  constant  coefficients  of  g  are  all  real  and  positive;  therefore,  the 
roots  of  g  =  0  are  all  complex.  This  can  be  seen  as  follows.  If  there 
were  another  real  root,  it  would  have  to  be  negative.  Let  y^  be  real  and 
posi tive.  For  y^  >  y^, 

g(-y^)  =  (y^^  +  y/yi)(y]^-yo)  *  .  (36c 

Similarly  for  > 

g(-yi)  =  +y^^  +  i>o.  (36d 

g  is  always  positive  for  a  real  y  and  g=0  can  have  only  complex  roots. 

Thus,  there  is  only  one  peak  in  the  reflectivity.  As  increases,  the 
peak  shifts  to  larger  5.  It  follows  that  if  the  semiconductor  is  saturated, 
the  optimum  place  for  phase  conjugation  is  slightly  above  the  band  edge, 
and  if  there  is  significant  absorption  at  .j,  the  place  to  be  is  slightly 
below  the  band  edge. 


ACKNOWLEDGEMENTS 


We  gratefully  acknowledge  the  support  of  AFSOR  for  this  work. 


2-24 


APPENDIX 


For  the  convenience  of  the  reader,  we  explicitly  display  the 
contour  integrations  which  yield  (7a)  and  (7b).  After  the  angular 
integrations  are  performed,  one  obtains  four  types  of  integrals  which 


can  be  written  as 


I^(A)  =  f  dy  y^(y^+A)‘^(y^+A*)'^  , 

oOO 

I2(A,B)  =  /  dy  y^(y^+A)‘^y^+B)‘^(y^+B*)’^  , 

.00 

12(A)  =  J  dy  y^(y^+A)‘^  , 

••00 

1^(A,B)  =j  dy  y^[(y^+A)^  +  B^]"’ 


The  integrand  of  can  be  written  as 


/^(y^+A)'^(y^+A*)‘^  =  A(A-A*)'^(y^+A)’^  -  A*(A-A*)‘^(y^+A*)'^ 


.  (A-A*)'^y^(y^+A)'^ 


An  integration  by  parts  gives 


f'^dy  y^(y^+A)-^  =  2'^/  dy(y^+A)'^ 


Thus , 


I,  (A)  -  2'^(A-A*)"^[(A+A*)  f  dy(y^+A)'^  -  2A*  |  dy(y^+A*)’'^] 

*  ^  _  m  —  QO 


(A-1) 


(A-2) 


(A-5) 


( A-6) 


(A-7) 
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Let  A  =  where  a  is  real  and  positive  and  -  ^  <  0  <  ^  .  There  are 

two  distinct  roots  for  +  A  *  0,  which  are  =  a  exp^i($+  1-)^  and  » 
a  exp|ji(<{)+  Closing  the  contour  in  the  upper  half, 

/■foo  - 

dy  (y^+A)"^  =  (ir/a)  exp  (-i^)  (A-8) 

.CO 

=  iT(l+A*/lAl)(2jAl+-A+A*}'^'^^  . 

Using  (A-8)  in  (A-7),  one  finds  that 

Ij(A)  =  (TT/4a^)(l+cos  2()))''^  exp  (A-9) 

=  (tt/2)(1+A*/1A1)(21A!+A+A*)‘^''^  . 

The  integrand  of  can  be  written  as 

y^(y^+A)'^(y^+B)"^(y^+8*)‘^  =  -A(A-B)'^(A-B*)’^(y^+A)'^  (A-10) 

-B(B-A)'HB-B*)‘^(y^+B)"^  -  B*(B*-A)'HB*-B)'^(y^+B*)'^  . 

Now  use  (A-8) ; 

l2(A,B)  =  -njA(l+A*/iAi){A-B)'^(A-B*)‘^(2iAi+A+A*)‘^''^  (A-11) 

+  8(l+B*/(B!)(8-A)'^{B-B*)'^(2lB]+B+B*)'^''^ 

+  B*(l+8/!Bj )(B*-A)‘^(B*-B)'^(2{Bl+B’»+6)"^^^  } 

=  ’T(A-B)'^(A-B*)'^[(A+)B5)(2lBi+B+B*)'^''^ 

-(A+iAi)(2|A(+A+A*)'^''^]  . 


2-26 


In  (A-3),  an  integration  by  parts  yields 


13(A)  =  2’Y  dy  (y^+A)-^  (A-U) 

^  CO 

=  (7r/2)(l+A*/|A|)(2|Al+A+A*)'^^^ 

Finally,  the  integrand  of  (A-4)  can  be  written  as 

y^[(y^-*-A)^+B^]'^  =  (i/2B)[(A-iB)(y^+A-iB)‘^  -  (A+iB) (y^+A+iB)’^]  ,  (A-13) 

where  A  and  B  are  real.  Now  using  (A-8)  yields 


14(4,8)  =  (Tr/8/2){AW)^/'^(/^^+B^  . 

Let  =  (47r^)'^(2m^/fi^)^^^.  Using  (A-7),  (A-11),  (A-12)  and  (A-14), 


one  finds 

|dl<(2Tr)'^  =  2^1^(liu;+E5+iY^y)  ,  (A-15a) 

|dl<(2TT)-2  L^^(u))!L^^(u))|2  *  2^I^(E5-tiui-iY^^)  , 

j'd^(ZTr)'^  L^^(ijo)  iL^y{aj)  =  ^0^2^ ’  (A-15c) 

j’dt(2TT)“^  ’  (A-15d) 

/di<(2Tr)'2[L^^(u))]2  =  EQl3(£g+1iw+iY^^)  ,  (A-15e)  | 

|d^(2T^)‘^[L^^(uj)]2  -  2Ql3(Eg-1l«-iY^^)  ,  (A-15f)  I 

J'd^(27T)-^lL^^(o.)l2  =  2^I^{Eg-1ia,.Y^^)  .  (A-159)  ' 

/di?(2TT)-^lL^^(i^)j^  =  2^I^(Eg+T)co.Y^y)  .  (A-15h) 

Collecting  all  of  these  together  gives  {7a)  and  (7b). 
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FIGURE  1. 
FIGURE  2. 

FIGURE  3. 

FIGURE  4. 


FIGURE  CAPTIONS 


Model  band  structure. 

The  geometry  of  four-wave  mixing,  and  A^  are  the 
amplitudes  of  the  conjugate  and  signal  waves, 
respectively.  A^and  A^  are  the  amplitudes  of  the 
pump  waves. 

=  oj+w-o))!  as  a  function  of  the  frequency 
detuning  per  linewidth. 

The  reflectivity  of  conjugate  waves  as  a  function  of 
the  frequency  detuning  for  various  optical  thicknesses. 
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ABSTRACT:  A  series  expansion  in  empty  resonator  modes  is  used  to  approximate  the 
field  In  a  cavity  filled  with  a  nonlinear  gain  medium.  The  time  dependence  of 
the  expansion  coefficients  may  be  obtained  using  Maxwell's  equations  and 
orthogonality  properties  of  the  modes.  The  equations  are  solved  numerical Iv  to 
determine  the  solutions  in  the  vicinity  of  the  stable  point.  The  effects  of 
changing  the  number  of  modes  and  the  acuracy  of  the  overlap  integral  are 
examined.  The  results  are  found  to  be  most  reliable  for  cavities  where  all  modes 
but  one  are  nearly  suppressed  and  stable  operation  is  possible.  Self-sustain’nq 
pulsations  between  transverse  modes  and  bistability  were  discovered  for  some 
conditions. 


I.  INTRODUCTION 


In  order  to  describe  the  time  dependent  properties  of  lasers  it  is 
natural  to  construct  a  set  of  first  order  equations  of  motion.  The  Maxwell 
equations  for  the  laser  fields  are  second  order  in  both  time  and  space 
derivatives  and  are  therefore  very  difficult  to  work  with  directly.  If  one 
assumes  that  the  fields  may  be  represented  by  a  time  dependent  expansion  in  free 
space  modes,  the  spatial  variations  of  the  fields  can  be  separated  by  insertinq 
this  expansion  into  the  nonlinear  Maxwell  equations.  Using  orthogonality 
properties  of  the  modes,  and  neglecting  second  order  effects,  there  results  a  set 
of  first  order  differential  equations  for  the  expansion  coefficients  of  the  laser 


The  free  space  modes  are  plane  waves,  which  does  not  lead  to  much  error 
for  lasers  with  stable  cavities  or  Fabry-Perot  resonators.  However,  high  energy 
systems  often  include  unstable  resonator  cavities^^*'^^  because  of  their  large 
mode  volumes,  excellent  discrimination,  and  efficient  outout  coupling.  The  modes 
of  these  cavities  are  quite  different  from  plane  waves,  and  may  be  obtained  by 
solving  the  Fresnel-Kirchoff  linear  integral  equation^^‘®\  For  such  cavities 
expansions  in  empty  resonator  modes  can  be  used  to  describe  the  spatial 
variations  of  the  fields,  and  time  dependent  behavior  may  be  treated  in  the  usual 
way  with  a  set  of  first  order  time  differential  equations  for  the  expansion 
coeff  icients^^\ 

The  careful  treatment  of  the  saturable  gain  problem  for  cw  operated 
unstable  resonators  involves  a  nonlinear  integral  equation  with  a  volume  rather 
than  a  surface  integral.  In  practice  this  approach  proves  to  be  intractible. 
Consequently,  the  equations  of  motion  resulting  from  the  use  of  expansions  in 
unstable  resonator  modes  should  be  useful  for  stable  operation  as  well  as  for 
time  dependent  properties.  In  this  oaoer  the  time  dependent  behavior  in  the 
vicinity  of  the  stable  point  is  examined. 
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In  Section  II  the  equations  of  motion  for  cylindrically  symmetric 
homogeneously  broadened  lasers  are  developed  with  explicit  expressions  for  the 
normalization  and  overlap  integrals.  The  resulting  equations  of  motion  are  in 
the  form  of  Lamb  Theory^^^,  and  are  in  a  sufficiently  general  form  to  include 
both  transverse  and  longitudinal  degrees  of  freedom.  The  numerical  evaluation  of 
the  overlap  integrals  is  outlined  in  Section  III.  Fast  Hankel  transforms  were 
used  to  generate  the  unstable  resonator  modes  in  the  interior  of  the  cavity,  and 
the  volume  integral  was  calculated  with  Simpson's  rule.  In  Section  IV,  the 
stable  solutions  to  the  equations  of  motion  are  described,  and  the  use  of  these 
soluti'  i  to  calculate  output  couplings,  frequency  pulling  effects,  power 
fluctuations,  and  line  broadening  is  discussed.  The  time  dependent  behavior  in 
the  vicinity  of  the  stable  point  is  treated  by  expressing  the  equations  of  motion 
in  this  region  in  the  simpler  Van  der  Pol  form.  In  Section  V,  numerical  results 
are  given.  An  M  =  2  cavity  was  used,  with  values  of  Ng^  at  both  a  peak  and  a 
crossing  point  of  the  iAj  vs.  Ngp  curve.  An  examination  of  the  accuracy  of  the 
overlap  integrals  and  the  necessary  number  of  modes  was  made.  It  was  found  that 
for  some  situations  the  solutions  of  the  equations  of  motion  did  not  stabilize, 
and  self-sustaining  pulsations  between  transverse  modes  appeared,  a  phenomenon 
which  cannot  be  explained  by  geometric  optics  descriptions.  When  one  mode  did 
achieve  dominance,  the  energy  distribution  of  the  cavity  was  little  affected  by 
nonlinear  gain,  although  the  frequency  pulling,  output  coupling,  and  linewidths 
were  changed  considerably  from  their  empty  resonator  values.  Finally,  we 
conclude  with  a  brief  discussion  of  the  use  of  exoansions  in  unstable  resonator 
modes  to  describe  nonlinear  systems. 
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II.  EXPANSION  IN  EMPTY  RESONATOR  MODES 


The  empty  resonator  modes  are  eigenvalue  solutions  to  the 
Fresnel-Kirchoff  integral  equation^^®"^^'.  The  kernel  to  this  linear, 
homogeneous  equation  is  non-hermitian  and  consequently  the  eigenvalue  solutions 
do  not  form  a  complete  set^^^^.  The  solutions  are  orthogonal and  may  be  used 
to  approximate  an  arbitrary  function  using  a  series  expansion  with  overlap 
integrals  to  compute  the  expansion  coefficients.  The  accuracy  of  this 
approximation  is,  in  general,  limited  due  to  the  non-completeness  of  the  empty 
resonator  modes.  However,  the  specific  function  we  wish  to  represent  is  the 
distribution  of  radiation  on  a  cavity  with  a  nonlinear  gain  medium  and  it  is  well 
known  that  this  distribution  is  not  very  much  different  from  the  lowest  loss 
empty  resonator  mode^^^"^®^.  Therefore,  one  might  expect  an  expansion  in  empty 
resonator  modes  to  give  reasonable  results  for  the  nonlinear  problem.  A 
relatively  small  number  of  terms  in  the  series  will  be  required,  by  contrast  with 
a  fourier  series  expansion  of  the  transverse  spatial  dependence,  which  has  a 
large  number  of  terms  due  to  the  effects  of  edge  diffraction.  Since  the 
difficulty  of  solution  increases  with  n^  where  n  is  the  number  of  terms  in  the 
series  approximation,  one  is  led  to  consider  the  use  of  empty  resonator  modes  to 
describe  nonlinear  gain  systems. 


The  electric  field  of  an  empty  unstable  resonator  may  be  expressed  as 
the  following  expansion  in  empty  resonator  modes: 


E(r,t)  =  Z  b„.  U„.  (r,z)e 


n£q 


n£q  n£q' 


(la) 


nP.q 


(r,z)e 


ikn£q(L/2-z) 
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n2q 


(lb) 


In  Eq.  (1)  ap  and  a0  are  forward  and  backward  amplitudes,  n  labels  the  radial,  ^ 
the  azimuthal  and  q  the  longitudinal  degrees  of  freedom.  The  cavity  is  asssumed 


to  be  cylindrical  with  equal  mirrors.  The  amplitudes  are  related  to  the  slowly 
varying  mirror  amplitude  by: 


X  Jj^(2TTNLxy/2)  , 
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ydye 
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ag"'^(x,z)  =  (-i)'"l  2.Nte 
X  J^(27'NLxy/2)  (y). 
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Here  x  and  y  are  radial  coordinates  normalized  to  the  mirror  radius,  M  is  the 
cavity  magnification,  N  is  the  q  dependent  cavity  Fresnel  number,  and  d'^^^'^fx)  is 
the  usual  solution  to  the  Fresnel-Kirchoff  integral  equation,  a  slowly  varying 
function  in  the  geometric  optics  limit.  Both  forward  and  backward  amplitudes 
have  been  explicitly  included  because  ap  ^  ag*  for  lossy  modes. 

The  k^pq  and  in  Eq.  (1)  are  complex  quantities,  determined  by  the 
eigenvalues  of  the  Fresnel-Kirchoff  integral  equation.  Each  mode  in  Eq.  fl) 
satisfies  Maxwell's  equations  in  the  axial  approximation  with  the  mirror  boundary 
conditions  Ein,if.^gp  =  0.  Since  has  a  positive  imaginary  part  for  all  lossy 
modes,  the  radiation  in  the  cavity  decays  with  time.  The  complex  part  of  kp.p 
compensates  for  the  decrease  in  the  magnitude  of  ap  as  a  wave  propagates  across 
the  cavity  in  the  forward  direction.  As  a  result  Eq.  H)  is  symmetric  under 
spatial  inversion.  If  some  radiation  is  introduced  into  the  cavity,  and  enough 
time  passes  to  form  the  mode  patterns,  the  initial  direction  of  the  radiation 
launched  into  the  cavity  is  forgotten.  The  radiation  distribution  is  spatially 
symnetric  about  the  midpoint  of  the  cavity,  and  exponentially  decreasing  in  time. 


To  account  for  time  dependent  properties  in  a  more  realistic  manner, 
the  expansion  coefficients  are  assumed  to  be  slowly  varying  functions  of 
time.  Since  b-.-  is  complex,  the  e^'^njiq^  term  can  be  absorbed  as  follows; 


^nsiq® 


Y  (t)e 


where  Yp()q  and  0p.q  are  real  functions  of  t  to  be  determined,  and  the 


instantaneous  frequency  of  the  mode  is  enjiq* 


The  resulting  expansion  is  substituted  into  the  wave  equation  for  the 


cavity, 


^2  1 

c  d  t 


e7  ^ 


where  P  is  a  complex  function  of  the  quasi-linear  form  P  =  g(ElE.  On  the  left 
side,  the  usual  axial  aporoximation^^^^  is  made  which  consists  of  neglecting 
second  derivatives  of  ap  3  with  respect  to  z.  Neglecting  also  second  derivatives 
with  respect  to  time,  and  taking  into  account  that  the  amplitudes  ap^5  satisfy 


Tz 


the  left  side  of  Eq.  (4)  is; 


)^f,b=o 


^^^neq  ■  niq*^  ^  '^'n£q  '  ^n£q  "^niq  j|'^nf,q 


;  qjb'  „  {r,2)e''^'^  ,  (6) 


where  Un,qfr,z)  is  the  spatial  part  of  the  mode,  given  by  Eq.  (Ibl, 
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Next,  two  further  approximations  are  made  which  are  consistent  with  the 
neglect  of  second  order  effects: 


(7a) 


•  • 


®n£q^niq  ^n£q 


(7b) 


Here  oj  is  a  real  frequency  equal  to  :Tqc/L,  where  q  is  the  longitudinal  mode 
number. 


The  right  side  of  Eq.  (4)  is  a  complicated  nonlinear  function,  and  for 
many  cases  the  linear  gain  is  large.  For  behavior  near  the  stable  point, 
saturation  is  a  factor,  and  y  g(Elc.  Then  the  right  side  of  Eq.  (4)  is  small 
to  first  order,  and  may  be  replaced  by 


g(E)£. 


(8) 


Next,  multiply  both  sides  of  the  resulting  equation  by 

n'i'q-  -ie'^i,'  '‘"n-rq-  -'’^n-y.'q- 


(9) 


and  integrate  over  all  taking  into  account  the  orthogonality  properties  of 
the  empty  resonator  modes.  Define  the  normalization  matrix  W  and  the 
polarization  matrix  P  by: 

ik  ...  ,(z-L/2) 
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where  <  >^  reoresents  an  average  over  the  cavity  volume.  The  matrices  W  and  P 
are  diagonal  in  8,  for  most  reasonable  functions  gfE). 


The  complex  form  for  k^j^qC  is 


"  “  *  Snilo  ♦  • 


(U) 


where  is  a  diffractive  frequency  shift,  Tj^^^q  the  diffractive  loss,  and  both 
B  and  Fc  are  «  w.  Finally,  the  cavity  equations  reduce  to  the  following  simple 
rate  equations: 
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Explicit  forms  for  the  normalization  matrix  W  and  t^e  polarization 
matrix  P  are  given  below.  The  orthogonal ity  relation  between  transverse  modes  of 
the  same  longitudinal  wave  is  given  by: 

nTq  ,  niq  ''^'^n'li’q  "  h  k  , 

<(ag  ap  e  e 


'^nn'^Zi' 


ydye 
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where  is  the  Fresnel-Kirchoff  solution. 


(13) 
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The  rapidly  oscillating  quadratic  term  in  the  integrand  results  from 
the  opposite  curvatures  of  the  forward  and  backward  waves.  In  order  to  evaluate 
the  interior  amplitudes  ap  and  ag  using  the  Fresnel-Kirchoff  solutions,  it  proves 
convenient  to  define  the  following  auxilliary  functions  in  Hankel  transform 
space: 


(Q)  ,  j'g-iirN(M-l)y2  ^n£qj  ^ 


(14a) 

(14b) 


where  =  0  outside  the  mirror,  for  x  >  1.  The  Hankel  transformed  interior 
amplitudes  are  given  by  Eq.  (22).  The  integral  in  Eq.  (13)  in  terms  of  these 
auxilliary  functions  -is: 


«nn' 
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Using  Eqs.  (2),  (10),  (14),  and  (15),  it  is  easy  to  derive  the  following  form  for 
the  W  matrix: 


Vi 
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where 


n'Jl'q' ,  n£q 
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(17b) 


For  large  q,  it  is  an  extremely  good  approximation  to  set  T  =  0  because  of  the 
rapidly  oscillating  z  dependence  of  the  integrand.  When  this  is  done,  the  theory 
is  identical  to  that  of  reference  (9).  Note  that  for  q’  =  q,  N'  =  N  and 
Eq.  (17a)  reduces  to  Eq.  (15)  due  to  the  orthogonality  prooerties  of  the  empty 
resonator  modes. 

The  evaluation  of  the  P  matrices  depends  on  the  specific  form  for  the 
nonlinear  susceptibility.  In  general,  the  calculation  involves  complicated 
overlap  integrals.  The  homogeneous  broadening  form  for  the  nonlinear 
susceptibility  is  used  here, 

g(E)  =  ig^  Cq  C/..0  (l+lE!^)'^  ,  (18) 


where  !E!^  has  been  normalized  to  the  saturation  intensity,  and  g^  is  the  linear 
gain.  Couplings  between  modes  of  different  i  are  neglected,  as  well  as  couplings 
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Jut. 


between  different  longitudinal  modes.  Interference  between  forward  and  backward 
waves  is  expected  to  be  significant  for  unstable  cavities^^®^.  The  frequency 
dependence  for  g(E)  is  not  specified,  since  the  laser  linewidth  is  assumed  to  be 
set  by  the  cavity  and  not  the  medium. 


As  in  Eq.  (16)  for  W,  the  integral  for  P  splits  into  two  terms. 
However  the  term  analogous  to  T  in  Eq.  (17b)  is  no  longer  negligible  because 
nonlinear  coupling  between  forward  and  backward  waves  leads  to  a  slowly  varying 
component  of  the  z-integral  and  hence  a  non-zero  contribution  to  the  integral. 
After  some  calculation,  the  effective  P  matrix  is  found  to  be 
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Here  the  T's  are  the  losses  as  in  Eq.  (11),  and  f  is  given  by: 
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Here  the  Indices  i,j  are  shorthand  for  the  mode  Indices  ntq.  Setting  B  =  0  in 
Eq.  (20)  is  equivalent  to  neglecting  the  effects  of  interference  between  forward 
and  backward  waves. 


1 


III.  EVALUATION  OF  OVERLAP  INTEGRAL 


The  solutions  to  the  Fresnel-Kirchoff  integral  equation  may  be  obtained 
by  well  known  methods^^"®^.  The  empty  resonator  modes  and  eigenvalues  used  here 
were  calculated  using  the  Butts-Avizonis  asymptotic  method,  with  the  following 
general  form^®^  substituted  for  their  F{x,t)  function: 


F(x,t)  =  K,[l]  =  f  j  (2tx)  , 

^  k=0 


where  t  =  -n-NM,  and  K[.]  is  the  Fresnel-Kirchoff  integral  transform.  The  mirror 
amplitude  c  is  given  in  terms  of  the  basis  functions  by 


♦  (X)  =  1  *  F„(x)  . 


Here  ap  is  obtained  by  finding  the  roots  of  a  N^^  -  order  polynomial  and 


Fp(x)  =  F(x/M" 


n  ou 

where  Mp  The  eigenvalues  and  eigenvectors  for  the  six  lowest  loss 
modes  for  both  the  Neq  =  6.5  and  Neq  =  7.02  M  =  2  resonator  were  obtained  and 
stored  on  disc  and  used  as  data  for  the  mode  mixing  calculation.  Only  the  ii  =  0 
case  was  treated.  The  losses  and  shifts  for  these  modes  are  obtainable  from  the 
eigenvalues  and  are  given  in  Table  I. 

Next  the  auxilliary  functions  in  Q-space  defined  by 
Eq.  (15)  are  calculated  using  the  Fast  Hankel  Transform  (FHT)  algorithm^^^^, 
which  approximates  the  Hankel  transform  by  a  finite  transform.  Since  the 
interior  values  of  the  amplitudes  are  determined  from  the  Fresnel-Kirchoff 
solution  by  the  application  of  a  finite  Hankel  transform,  the  FHT  is  particularly 
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appropriate  for  this  application.  The  conditions  for  the  accuracy  of  the  FHT  are 
well  known^^^),  and  lead  to  the  following  mesh  for  x  and  Q  space: 


"n  =  "o® 


an 


Q  =  Q  e 


an 


a  =  1/4S 
-aH 

"o  -  ® 


(25) 


where  6  is  the  highest  spatial  frequency  component  considered,  chosen  here  to  be 
2.5  lowest  truncation  point  is  less  than  1/4  of  the  highest  spatial 
frequency,  and  M  =  512  is  the  number  of  points  in  the  mesh. 


Once  Op(Q)  and  413(0)  are  known,  the  Hankel  transformed  interior 
amplitudes  are  simply: 


ap^B(q,2)  =  e 


n'TtO  2/NL 
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(26) 


The  transverse  spatial  dependence  of  the  interior  amplitudes  is  found  by  again 
Hankel  transforming  Eq.  (26).  For  each  interior  point  z,  the  amplitude  is  given 
as  an  M-vector.  Finally  f{x,z)  is  determined  from  Eq.  (20)  and  the  double 
integral  in  Eq.  (19)  is  evaluated  using  the  trapezoidal  rule.  It  is  only 
necessary  to  integrate  halfway  across  the  cavity  due  to  the  symmetry  of  the 
problem. 


The  number  of  points  in  the  z  axis  grid  does  not  have  to  be  very  large. 
Numerical  solutions  of  the  nonlinear  integral  equation  assume  that  the  nonlinear 
gain  is  concentrated  on  a  small  number  of  thin  slabs,  and  free  space  propagation 


is  used  in  the  intervals  between  the  slabs.  Free  space  propagation  is  given  bv 
Eq.  (26),  and  it  may  easily  be  seen  that  the  use  of  the  trapezoidal  rule  with 
Osiab  points  in  the  z  grid  is  just  as  accurate  as  the  evaluation  of  the  nonlinear 
Fresnel-Kirchoff  integral  transform  using  the  same  number  of  thin  slabs.  The 
single  thin  slab  approximation  has  been  used  extensively  for  low  loss 
cavities^^^"^^\  and  has  gives  reasonable  results  in  this  regime.  Consequently, 
a  relatively  small  number  of  gain  slabs  should  be  sufficient  to  describe  the 
radiation  in  more  lossy  cavities.  For  the  present  calculation  we  used 
5  <  njigjj  <  45,  and  making  the  grid  finer  did  not  change  the  results 
significantly,  (see  Figure  5). 


IV.  SOLUTIONS 


CW  laser  operation  may  be  adequately  described  by  the  stable  solution 
to  the  equations  of  motion.  The  time  dependence  of  the  fields  in  the  immediate 
vicinity  of  the  stable  point  specifies  such  quantities  as  the  frequency 
fluctuation,  power  fluctuation,  and  the  field  spectrum.  Consequently,  we  first 
consider  the  stable  solutions  to  Eq.  (12). 

The  determination  of  the  stable  point  involves  solving  a  set  of 
coupled,  nonlinear  equations.  This  was  done  in  the  present  case  by  setting 
v^-  =  0,  rewriting  the  equations  in  the  form  x  =  ffxl,  then  using  iteration, 
Xj+1  “  This  nonlinear  mapping  can  lead  to  limit  cycles  or  chaotic  behvior 
rather  than  a  stable  point  in  some  cases^^^^.  The  stable  solutions  to  the 
resulting  n  complex  equations  give  the  n  T^'s,  the  n-1  relative  phases  -  “-j.i, 
and  the  transverse  frequency  shift  Aw,  for  a  total  of  2n  real  variables.  Only 
the  relative  phases  may  be  determined  because  the  nonlinear  term  contains  E  in 
the  form  !e|^,  which  is  invariant  to  overall  phase  translations.  The  frequency 
shift  Aw  is  the  lockel  value  of  the  transverse  frequencies,  and  is  related  to  the 
phases  by  Si  -  ^  -  Aw. 

It  may  seem  surprising  that  the  frequencies  should  lock  together,  since 
the  transverse  frequencies  are  often  spaced  further  apart  than  the  longitudinal 
frequency  separation  c/2L,  and  the  susceptibility  is  purely  imaginary,  at  each 
point  in  space.  However,  an  effective  nonlinear  index  of  refraction  arises  from 
the  coupling  of  edge  diffracted  rays  in  the  center  of  the  cavity.  Due  to  the 
complicated  phase  dependence  of  the  individual  modes,  the  overlap  integrals  in 
Eq.  (12)  are  not  purely  imaginary,  even  though  the  susceptibility  is.  The 
overlap  integral  coupling  between  modes  is  non-local  in  nature,  and  this  type  oT 
coupling  can  lead  to  frequency-pulling  effects  even  if  the  index  of  refraction  is 
locally  linear.  Since  changing  the  relative  phases  changes  :En  and  therefore 
affects  the  overlap  integrals,  transverse  mode  locking  can  occur. 


This  type  of  mode  locking  is  different  from  the  case  where  longitudinal 
mode  frequency  differences  lock  together  leading  to  the  formation  of  pulses^^^. 
For  unstable  cavities,  the  transverse  frequencies  are  spread  so  irregularly  that 
it  is  difficult  for  population  pulsations  at  the  frequency  differences  to  appear. 
If  there  is  coupling  between  a  number  of  different  longitudinal  modes,  a  case  can 
be  envisioned  where  all  the  transverse  modes  of  each  longitudinal  mode  lock 
together,  and  the  differences  between  longitudinal  frequencies  then  lock  together 
by  the  usual  means. 

Once  the  stable  solution  for  and  is  obtained,  the  mirror 

amplitude  is  given  by  <i)(x)  =  E  Y,e’®i<i)i(x)  and  the  locked  frequency  of  the 

i  '  ' 

transverse  modes  by  w  +  iw.  The  output  coupling  for  the  stable  solution  may  be 
determined  from  the  effective  eigenvalue,  defined  by 

,  -  ^  >  (27) 

eff  ■ 

where  t  is  the  converged  mirror  amplitude,  the  brackets  represent  a  mirror 
average,  and  K  is  the  linear  Fresnel-Kirchoff  integral  transform  which  propagates 
the  solution  once  across  the  empty  cavity.  Using  the  expansion  of  ?  in  empty 
resonator  modes  and  orthogonality  properties,  one  finds 


22ie.  „ 

'<-i  > 

1  ' 


(23) 


where  are  the  eigenvalues  of  the  empty  resonator  modes.  The  losses  are  given 
by  1  - 

To  study  the  fluctuations  about  the  stable  ooint,  a  set  of  Langevin 
equations  may  be  formed  by  adding  random  noise  terms  on  the  right  side  of 
Eq.  (12).  Alternatively,  Eq.  (12)  may  be  used  to  wr’te  down  the  Fokker-Plank 
equation  for  the  system(^^\  Due  to  the  complicated  nature  of  the  right  hand 
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side  of  Eq.  fl2),  neither  of  these  approaches  is  useful  unless  one  is  able  to 
solve  the  resulting  equations  numerically. 


Here  we  take  another  approach  which  is  approximate,  yet  provides 
qualitative  insight  into  the  behavior  of  the  system  as  well  as  yielding  useful 
numerical  predictions.  The  equations  of  motion,  Eq.  fl2),  in  the  vicinity  of  the 
stable  point  are  approximated  by  a  simpler  form  which  may  be  treated 
analytically,  that  of  a  Van  der  Pol  oscillator^-': 

■^i  =  ‘"i^i  -  ?i.i^i'^i-  (29a) 


-  CO  = 


(29b) 


Here  and  are  synmetric,  and  the  relative  phases  are  assumed  to  be  fixed.  The 
3^-j,  and  parameters  were  chosen  by  numerically  evaluating  and 
using  Eq.  (12)  for  a  large  number  of  randomly  chosen  in  the  immediate 
vicinity  of  the  stable  point,  and  then  least  squares  fitting  to  the  form  Eq.  129^ 
with  the  restrictions  that  Yi  =  0  and  9,-  =  *  lu  at  the  stable  point. 

The  Fokker-Plank  method  is  used  here,  and  the  noise  terms  added  on  the 
right  hand  side  of  Eq.  (29)  are  assumed  to  be  Markovian: 


<n./t) 

n,^(t' 

)  >-  D  i(t-t’ ) 

(30a) 

<ng!t) 

)>=  R  3(t-f)  . 

(30b) 

This  assumotion  leads  to  flat  spectra  for  the  power  fluctuations,  freouencv 
fluctuations,  and  field.  To  obtain  more  realistic  oredictions  for  the  soect^al 
behavior  or  these  quantities,  a  Langevin  equation  treatment  should  be  used. 
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V 


The  next  approximation  is  to  consider  only  the  amplitude  noise  term, 
Eq.  (30a).  By  neglecting  phase  fluctuations,  the  Fokker-Plank  equation  can  be 
solved  trivially  since  Eq.  f29a)  may  be  derived  from  a  potential.  Amplitude 
fluctuations  will  lead  to  frequency  shifts  through  Eq.  (29b),  which  is  assumed  to 
react  instantaneously  to  small  amplitude  fluctuations. 

With  these  assumptions,  the  solution  to  the  Fokker-Plank  equation  is 

P(Yi...Yn)  =  exp  |(2aY^  -  Y^By^)/2d|,  (31) 

where  y  is  the  vector  y^  ,  is  the  vector  a^-,  6  is  the  matrix  6^-j,  and  N  is  a 
normalization  constant.  Eq.  (31)  can  be  used  to  obtain  estimates  for  the  power 
fluctuations  by  assuming  that  all  of  the  amplitudes  except  the  largest  is  fixed 
at  its  stable  value,  and  finding  the  point  where  P  drops  1/2  of  its  value  at  the 
stable  point.  The  FWHM  of  the  power  fluctuations  is  approximately: 


AI/I  %  2Ay^/y^ 


8Dln2 

Y  ®ii 


(32) 


where  yi  is  the  largest  mode  amplitude. 


To  obtain  a  rough  estimate  for  the  frequency  broadening  due  to 
amplitude  fluctuations,  only  those  fluctuations  which  maintain  the  frequency 
locking  of  the  modes  are  assumed  to  be  relevent.  The  locked  frequency  itself  is 
buffeted  about  by  these  fluctuations.  This  assumption  enormously  simplifies  the 
calculations  and  is  expected  to  give  reasonable  results  for  the  frequency 
broadening.  In  fact,  results  for  the  field  spectrum  E((3j)  obtained  using  this 
method  agreed  well  with  a  more  exact  calculation  made  for  the  cases  of  1,  2,  and 
3  modes.  Requiring  -  w  to  be  a  single  constant  determines  a  point  in  Y^  space 
through  Eq.  (29b),  and  the  multidimenional  distribution  P  in  Eq.  (31)  is  then 
reduced  to  a  one-dimensional  dependence  on  e^-  -  oj  =  6a): 


P(( 


■1 


-1 


=  N  exp  I (-Zap  u  -  uc 


■^3d‘^u)/2D 


(33) 


3-20 


Here  p  and  s  a'"e  the  matrices  and  a  is  the  vector  a^,  and  u  is  the 

vector  §,•  -  0)  -  =  5oj  -  a^.  The  broadening  is  obtained  by  determining 

the  full  width  at  half  maximum  of  P(icj),  and  is  easily  seen  to  be: 


Aw  = 


I  8D1n2 

V  Z  (p'^ep*^).. 

1*  1  '  J 


(34) 


It  is  necessary  to  give  some  value  for  D  in  order  to  obtain  numerical 
results.  The  noise  term  arises  from  spontaneous  emissions,  so  one  expects  D  to 

be  small  for  large  coherent  intensities  and  large  if  the  lifetime  of  the  upper 

laser  state  is  small.  Also,  0  should  be  directly  proportional  to  the  number  of 
atoms  in  the  upper  laser  state.  To  determine  D,  one  calculates  the  total  energy 

stored  in  a  Van  der  Pol  oscillator  driven  by  a  random  noise  term,  nYftl^^'^^.  The 

total  energy  is  obtained  in  terms  of  the  power  spectrum  of  n^ft),  Wj^fw).  When 
this  energy  is  set  equal  to  the  energy  in  the  cavity  due  to  spontaneous  emission, 
W_  is  determined  and  D  is  then  obtained  from 

ny 

D  =  TrW^^(w^)  ,  (25) 


where  Wj  is  the  laser  frequency.  The  result,  taking  into  account  the  particular 
definitions  of  the  y-jtj  and  the  fact  that  Eq.  (29)  was  obtained  from  a 
calculation  which  used  the  dimensionless  time  x  =  ct/2L,  is: 


D  = 


— 2 - 

4TT‘^t2LMl5 


(36) 
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Here  D  is  dimensionless,  as  are  a,  B,  and  x,  N2  is  the  number  density  of  the 
upper  laser  state,  t2  is  the  lifetime  of  the  upper  laser  state,  L  is  the  cavity 
length,  N  the  cavity  Fresnel  number,  and  is  the  saturation  intensity.  With 
the  above  relation  for  D,  the  power  fluctuation  and  frequency  broadening  may  be 
estimated. 


V.  RESULTS 


The  main  purpose  of  the  calculations  was  to  examine  the  usefulness  of 
the  idea  of  expressing  the  field  of  a  loaded  unstable  resonator  as  a 
superposition  of  empty  resonator  modes.  In  all  cases,  coupling  between 
transverse  modes  only  was  considered.  Since  only  one  longitudinal  mode  was 
present,  this  corresponds  to  what  is  usually  referred  to  as  single  mode 
operation.  By  gradually  increasing  the  number  of  transverse  modes  until  the 
resulting  spatial  distributions  and  frequencies  of  the  stable  solutions 
converged,  it  was  found  that  a  relatively  small  number  of  modes  were  needed  for 
an  accurate  description.  In  addition  both  the  transverse  and  the  longitudinal 
accuracy  of  the  overlap  integral  was  varied,  with  grid  density  increasing  until 
the  results  converged.  Finally,  an  examination  of  the  behavior  in  the  vicinity 
of  the  stable  point  was  examined  by  using  fits  to  Eq.  (291.  From  these,  the 
DOwer  and  frequency  fluctuations  were  estimated. 

Two  separate  sets  of  calculations  were  made  for  cavities  with  M  =  2, 
with  Neq  =  6.5  or  7.02.  On  the  plot  of  eigenvalues  versus  equivalent  Fresnel 
number  Neq  for  an  M  =  2  resonator,  there  is  a  peak  at  Neq  *  6.5  and  a  minimum  at 
Neq  =  7.02.  At  the  peaks,  which  occur  near  half  integral  values  of  Neq,  the 
lowest  loss  mode  is  well  separated  from  the  other  modes.  By  contrast,  the 
minimums  occur  near  integral  values  of  Neq  and  result  from  the  crossing  of  two 
separate  empty  resonator  modes  having  the  same  loss.  The  two  sets  of  conditions 
should  pretty  well  characterize  the  response  of  an  M  =  2  cavity  to  a  nonlinear 
gain  medium.  The  magnitudes  and  phases  of  the  eigenvalues  for  the  Neq  =  6.5  and 
Neq  =  7.02  empty  resonator  modes  are  given  in  Table  I. 

First,  the  dependence  of  the  stable  solution  on  gain  was  examined.  The 
results  for  the  Neq  =  6.5  three  mode  cavity,  with  a  very  coarse  integration 
mesh  =  5  are  shown  on  Figure  1  in  the  form  of  a  plot  of  the  area-averaged 
mirror  intensity  vs.  gain.  Three  modes  were  included  because  the  second  and 
third  states  have  equal  losses  and  thus  should  be  considered  together.  As 
expected,  the  lowest  loss  mode  has  the  largest  coefficient,  although  as  gain 
increases  the  coefficients  of  the  lossy  modes  become  somewhat  larger.  Also  shown 
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on  Figure  1  are  the  geometric  optics  area- averaged  mirror  intensities  obtained  by 
solving  the  Rigrod  equations  for  the  cavity,  with  couplings  between  forward  and 
backward  waves  included^^^*^^^ .  Note  that  the  threshold  for  the  geometric  optics 
case  occurs  at  g^L  =  1.39,  corresponding  to  the  geometric  optics  eigenvalue 
]X|  3  a  .5.  The  threshold  for  the  three-mode  solution  with  full  transverse 
field  variation  occured  at  the  lower  value  ggL  *  .93,  corresponding  to  the 
eigenvalue  of  the  lowest  loss  mode  U!  =  .627.  All  three  modes  are  less  lossy 
than  the  geometric  optics  predictions  because  diffraction  couples  radiation  from 
the  edges  back  into  the  center  of  the  cavity. 

The  geometric  optics  calculation  results  in  a  larger  average  intensity 
than  is  found  in  the  more  careful  three-mode  treatment.  One  reason  for  this  is 
that  as  gain  increases,  the  lossy  modes  become  excited  and  therefore  output 
coupling  increases.  On  the  other  hand,  output  coupling  remains  independent  of 
gain  for  geometric  optics,  and  therefore  one  expects  larger  fields  in  thi'  case. 
Also,  the  area  averaging  weights  the  edge  region  disproportionately,  where  the 
geometric  optics  wave  is  markedly  larger  than  the  lossy  modes.  In  any  case,  the 
geometric  optics  results  give  qualitatively  reasc.iable  predictions,  as  long  as 
there  actually  is  a  stable  solution  to  the  problem. 

As  gain  is  increased  further,  oscillations  eventually  appear  in  the 
three-mode  solution  and  the  iteration  does  not  converge,  even  though  the 
geometric  optics  problem  has  stable  solutions  for  arbitrary  gain.  This  behavior 
is  much  more  obvious  for  the  Neq  »  7.02  cavity,  where  there  are  two  modes  of 
equal  loss  competing  for  gain.  For  this  cavity  a  large  number  of  calculations 
were  made  for  the  two  mode  problem  with  n^i^ij  =  10  at  different  levels  of  gain. 
Since  these  two  modes  have  equal  losses,  the  problem  cannot  be  further 
simplified.  To  check  that  two  modes  were  enough,  calculations  were  done  for 
various  gains  with  3,  4,  and  6  modes,  and  the  results  were  found  to  be 
qualitatively  the  same  as  for  the  two  mode  problem,  because  all  the  higher  order 
modes  had  small  coefficients. 
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At  least  one  stable  solution  was  found  for  each  linear  gain.  The 
converged  values  of  and  Y2  for  g^L  <  8  are  given  in  Figure  2.  Near  threshold, 
Yj  is  largest.  As  gain  is  increased,  Y2  increases  continuously,  and  near  g^L  =  6 
Y2  suddenly  becomes  the  dominant  mode.  Above  this  point,  and  at  least  until 
g^jL  =  150,  Y  j  is  nearly  suppressed.  In  the  region  5.7  <  g^L  ^  6  two  separate 
stable  solutions  were  found,  implying  the  existence  of  bistable  operation  for 
this  cavity.  It  is  interesting  to  note  that  for  unstable  cavities,  a  locally 
nonlinear  index  of  refraction^^^^  is  not  needed  in  order  for  bistability  to 
occur.  The  discontinuity  in  the  solutions  near  g^L  =  6  was  not  found  for  the 
Ngp  =  6.5  cavity  and  is  typical  of  bistable  systems.  The  discontinuity  is  also 
evident  in  the  mirror  averaged  intensities  as  gain  is  increased  or  decreased 
through  the  bistable  regime.  It  is  expected  that  there  is  a  regime  of  stationary 
solutions  which  joins  the  two  branches  of  stable  solutions,  but  are  unstable  to 
small  perturbations.  Our  algorithm  was  not  designed  to  obtain  these  types  of 
solutions. 


In  addition  to  stable  solutions,  self-sustaining  pulsations  were 
discovered  in  the  region  of  bistability.  For  some  initial  conditions,  the 
interative  procedure  did  not  converge  to  a  stable  point  but  to  a  limit  cycle. 
Limit  cycles  for  g^L  =  6  and  g^L  *  7  are  displayed  in  Figures  3  and  4.  Similar 
limit  cycles  were  also  found  for  the  four  mode  problem,  so  the  phenomenon  is  not 
a  result  of  the  two  mode  restriction.  The  existence  of  a  limit  cycle  imolies 
that  periodic  self-sustaining  pulsations  between  the  two  lowest  loss  modes  occurs 
as  the  solution  traverses  the  limit  cycle.  It  is  interesting  that  either  a 
stable  solution  or  pulsating  behavior  can  occur  for  the  same  cavity  at  a  specific 
gain,  depending  on  the  initial  conditions. 

As  gain  is  increased,  the  effective  coupling  between  modes  increases. 
For  small  gains,  both  modes  can  oscillate  at  the  same  time,  but  as  coupling  is 
increased,  a  point  is  reached  where  one  mode  is  suddenly  suppressed.  This 
evidently  occurs  near  g^L  =  6.  For  larger  gains,  the  solution  with  both  modes 
oscillating  is  not  stable  and  this  is  precisely  where  we  find  the  limit  cycles 
representing  periodic  pulsations  between  the  two  modes.  Both  modes  are  trying  to 
oscillate  at  once,  and  the  strong  coupling  leads  to  instability. 
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Next,  the  integration  accuracy  and  the  effects  of  including  more  modes 
into  the  calculation  were  examined.  Figure  5  plots  the  mirror  amplitude  (x1 
for  a  cavity  with  Neq  =  6.5  with  the  number  of  longitudinal  integration  points 
set  at  =  5,  15,  25,  and  45.  Also  given  in  Table  II  are  the  amplitudes  and 

phases  for  various  integration  meshes.  Evidently,  increasing  the  accuracy  from 
’’slab  ■  leads  to  little  noticeable  change  in  the  results.  As  expected  from 

previous  thin  slab  calculations  where  as  little  as  one  slab  was  used,  great 
accuracy  of  integration  in  the  longitudinal  direction  is  not  needed.  For  all  of 
these  calculations  512  transverse  points  were  used.  These  are  not  spread 
uniformly,  but  are  concentrated  near  the  center  of  the  cavity  due  to  the 
exponential  Gardiner  transformation  used  in  the  FHT  algorithm^^®^ .  This  number 
of  points  was  adequate  to  treat  the  transverse  field  dependence  of  cavities  with 
N  10, 


It  was  expected  that  a  small  number  of  modes  would  be  necessary  to 
describe  the  nonlinear  behavior,  and  indeed  for  most  cases  three  to  four  modes 
gave  essentially  the  same  results  as  six  modes.  Table  III  gives  the  amplitudes 
and  phases  for  both  cavities  considered  at  gQl  =  5  with  the  number  of  modes 
increased  from  1  to  6.  The  modes  are  ordered  according  to  loss,  and  it  is 
evident  that  the  more  lossy  modes  have  smaller  amplitudes.  For  the  Neq  =  6.5 
cavity,  the  lowest  loss  mode  is  by  far  the  largest  with  about  94%  of  the  energy. 
The  >5  parameter  for  Neq  =  6.5  appears  to  be  quite  large,  however  this  represents 
only  about  3%  of  the  energy.  For  the  Neq  =  7.02  cavity,  the  first  mode  has  59% 
and  the  second  36%  of  the  energy,  more  evenly  distributed  since  the  coupling  is 
not  quite  strong  enough  at  g^L  =  5  to  suppress  either  of  the  lowest  loss  modes. 

To  demonstrate  the  effect  of  gain  on  the  distribution  of  the  cavity 
radiation.  Figures  6  and  7  compare  the  six  mode  g^L  =  5  result  for  i?  with  the 
empty  resonator  g^L  =  0  results.  The  radiation  in  the  Neq  =  6.5  cavity,  with  94% 
of  the  energy  in  the  lowest  loss  mode,  is  little  affected  by  the  gain.  The  main 
difference  between  the  curves  is  that  the  distribution  in  the  saturable  gain 
medium  is  more  strongly  suppressed  near  the  mirror  edge.  The  bumps  and  valleys 
of  the  ci>’"ves  match  quite  well,  showing  indeed  that  the  empty  resonator  modes 
give  a  good  description  of  loaded  cavity  behavior.  On  the  other  hand,  the 


labs  were 


-110.50 


The  six  mode  mirror  amplitude  ii)(x)|  for  the  =  7.02  g^L  =  5 
cavity  with  45  points  in  the  z-integration  grid  is  compared  with 
one  of  the  lowest  loss  modes  of  the  corresponding  empty  resonator 
(dotted  line). 


Neq  «  7.02  solution  exhibits  significant  deviations  from  the  empty  resonator  mode 
as  expected  since  the  two  lowest  loss  modes  are  of  equal  importance.  Evidently 
the  empty  resonator  modes  best  characterize  loaded  cavities  when  a  single  mode 
absorbs  most  of  the  energy. 

Next,  the  effect  of  nonlinear  gain  on  frequency  pulling  and  output 
couplings  was  considered.  These  quantities  are  determined  by  cavity  geometry  and 
diffractive  effects  for  empty  resonators.  The  frequency  shift  of  the  Neq  =  6.5 
cavity  with  g^L  =  5  was  -.20  c/2L,  and  the  Neq  =  7.02  cavity  was  frequency 
shifted  a  much  larger  value  =  -.51c/2L.  In  the  first  case  the  frequency 
pulling  is  much  larger  than  the  empty  resonator  diffractive  frequency  shift  found 
in  Table  I.  The  large  frequency  shift  for  the  Neq  =  7.02  case,  equivalent  to  1/2 
of  the  longitudinal  spacing,  is  typical  for  cavities  operating  near  a  crossing  of 
eigenvalues.  It  is  apparent  that  the  saturable  gain  medium  can  affect  frequency 
pulling  at  least  as  much  as  diffraction. 

The  effective  output  couplings  for  the  two  cases,  given  by  Eq.  (24), 
we'"e  72.8%  for  the  Neq  =  6.5  cavity,  and  66.3%  for  the  Neq  *  7.02  cavity.  By 
comparison,  the  geometric  optics  loss  in  both  cases  is  75%.  Note  that  the 
Neq  =  7.02  cavity  has  losses  determined  almost  entirely  by  the  two  crossed  modes 
with  64%  losses.  On  the  other  hand,  the  Neq  =  6.5  lowest  loss  is  only  61%,  while 
the  effective  loaded  cavity  losses  are  73%.  Couoling  to  the  higher  order  modes 
greatly  increases  the  losses  in  this  case,  and  it  is  not  understood  why  a  similar 
effect  does  not  occur  for  the  two  mode  Ng^  =7.02  cavity. 

To  determine  the  fits  to  the  Van  der  Pol  form  in  the  vicinity  of  the 
stable  points,  ^^-'s  were  varied  randomly  in  a  region  within  1%  of  their  stable 
values,  and  the  phases  were  kept  fixed.  For  each  set  of  amplitudes  chosen 
this  way,  the  and  were  calculated  using  the  overlap  integral  procedure.  A 
large  number  of  these  points  enabled  the  B,  o,  d  parameters  to  be  determined 
by  least  squares  fitting. 


1 
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For  this  procedure  to  give  reasonable  results,  the  amplitudes  of  the 
modes  should  vary  less  than  1%  from  their  stable  point  values  when  noise 
fluctuations  are  allowed  for.  From  Eq.  (32)  it  can  be  seen  that  in  order  for  the 
amplitude  fluctuations  to  be  less  than  1%,  0  must  be  less  than  10^.  D  is 
determined  from  laser  parameters  by  Eq.  (36),  and  using  some  typical  parameters 


y  =  10.6 
Nj  =  6.0  Torr 
t2  =  2.2  X  10"^  sec 
N  =  8.67 
1 5  =  1  watt  -  cm'^ 

L  =  100  cm 


(37) 


results  in  the  value  0  =  1.34  x  10"^,  which  implies  fluctuations  much  less  than 

1%. 

Consider  first  the  case  of  only  two  competing  modes.  The  fitted 
parameters  for  the  two  cases  at  g^L  =  5  are  given  in  Table  6.  The  stability 
analysis  of  two  mode  operation  involves  the  effective  gain  parameters  and 
ct2',  and  the  coupling  constant  C  defined  by^^^: 

”  ^2  ”  ^21*^1'^^11  (38) 

C  =  ^i2^21'^^11^22 


For  the  Neq  =  6.5  cavity,  the  parameters  in  Eq.  (38)  are  given  by 
C  =  2.23  X  10"'^,  =  .307,  and  =  .212.  Since  C<1,  linear  stability 
analysis  implies  the  existence  of  a  stable  point  with  both  modes  oscillating  at 
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Yi  »  1.526  and  Yp  “  .170.  This  Is  the  stable  point  that  is  actually  found.  On 
the  other  hand,  the  Ng^  =  7.02  cavity  has  C  =  4.43,  »  -.520,  and  -  .249. 
This  is  strong  coupling,  and  the  two  possible  stable  solutions  each  have  one  mode 
totally  suppressed.  The  actual  solution  found  has  both  modes  oscillating  with 
Yj  =  .832  and  2  ~  *^71.  According  to  linear  stability  analysis,  this  solution, 
though  stationary,  is  not  stable  with  respect  to  small  perturbations.  The  fact 
that  our  solution  is  Indeed  stable  indicates  that  the  fit  to  the  Van  der  Pol  form 
has  limited  usefulness  when  both  modes  are  oscillating.  For  larger  gains,  one  of 
the  Ngp  =  7.02  modes  was  indeed  suppressed  and  the  Van  der  Pol  fits  were 
consistent. 


The  Neq  =  6.5  cavity  with  g^L  =  5  appears  to  be  well  described  by  the 
Van  der  Pol  fit  in  the  vicinity  of  the  stable  point.  The  parameters  for  the  1, 
2,  3,  and  4  mode  calculations  for  this  case  are  given  in  Table  6.  With  D  chosen 
as  D  =  10"^  power  fluctuations  are  found  to  be  of  magnitude  AI/I  'v  10"'*. 

The  linewidths  are  given  by  the  F.W.H.M.  values  of  the  Pf6w)  curves, 
and  may  be  found  directly  from  Eq.  (31).  The  results  for  both  the  linewidths  and 
frequency  pulling  are  given  in  Table  V.  Increasing  the  number  of  modes 
considerably  changes  the  single  mode  results,  even  though  over  90iK  of  the  energy 
is  concentrated  in  the  lowest  loss  mode.  The  frequency  broadening  effects  do  not 
change  much  as  the  number  of  modes  becomes  greater  than  1,  indicating  that  a 
description  with  just  a  few  modes  gives  reasonably  consistent  results.  Note  that 
The  frequency  broadening  effects  considered  here  arise  from  amplitude 
fluctuations  and  mode  couplings,  not  phase  fluctuations. 
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Table  V.  Frequency  pulling  and  frequency  broadening,  in  units  of 
c/2L,  for  the  =  6.50  cavity.  The  shift  is  detennined 

by  the  stable  solution,  and  the  broadening  by  the  time 
dependent  behavior  in  the  vicinity  of  the  stable  point  under 
the  influence  of  a  random  noise  term. 


n  = 

1 

n  = 

2 

n  = 

3 

n  = 

4 

Pul  1 ing  Aw(c/2L) 

Broadening  Aa)(c/2L) 

-.072 

4.39  X  10‘® 

-.022 

33.24  X  10*^ 

-.023 

20.32  X  10'® 

-.018 
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o 

X 
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»— • 
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VI.  CONCLUSIONS 


The  use  of  series  expansions  in  empty  resonator  modes  to  describe 
lasers  with  nonlinear  gain  media  has  been  examined.  The  method  hinges  on  the 
evaluation  of  the  overlap  integrals,  which  involve  double  integrals  over 
longitudinal  and  transverse  dimensions.  It  has  been  found  that  ten  or  twenty 
grid  points  give  reasonable  results  for  the  longitudinal  integration,  while  a 
much  larger  number  of  points  are  needed  for  the  transverse  integrations.  This  is 
in  accordance  with  thin  slab  approximations,  which  have  been  shown  to  give  good 
results  for  a  relatively  small  number  of  slabs. 

Even  though  the  overlap  integrals  are  difficult,  this  method  is  still 
faster  than  direct  solution  of  the  nonlinear  integral  equation.  Assuming  that 
the  empty  resonator  modes  are  known  in  advance,  the  overlap  integrals  for  one 
iteration  involve  x  x  Ngigi,  complex  additions.  By  contrast, 

one  iteration  of  the  nonlinear  integral  equation  involves  3  x  N^-jgij  fast  hankel 

transforms,  or  6  x  Njigjj  x  ^transverse  ^  ^®92(^^transverse^  complex 
multiplications,  and  an  equivalent  number  of  complex  additions.  Since 
multiplications  are  about  ten  times  more  time  consuming  than  additions  and 
"modes  ^  evident  that  the  series  expansion  approach  using  overlap 

integrals  can  be  much  more  numerically  efficient  than  direct  solution  of  the 
nonlinear  integral  equation.  Provided  the  series  approach  also  gives  results 
equivalent  to  the  more  exact  integral  equation  method,  a  significant  increase  in 
efficiency  may  be  achieved  together  with  more  direct  physical  interpretation  of 
the  solutions. 

The  results  appear  to  be  best  for  cavities  operating  far  from  the 
crossing  points  of  the  eigenvalue  plot.  In  this  regime,  the  lowest  loss  mode  is 
well  separated  and  the  series  expansion  gives  reasonable  results  with  the  largest 
part  of  the  energy  is  concentrated  in  only  one  mode.  A  three  or  four  mode 
description  was  sufficient  to  characterize  the  response  of  the  resonator  to  a 
nonlinear  gain  medium,  and  the  radiation  pattern  was  qualitatively  and 
quantitatively  quite  similar  to  that  of  the  empty  resonator  mode.  By  contrast, 
the  effective  output  coupling,  frequency  pulling,  and  linewidth  varied 
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considerably  from  the  empty  resonator  results.  The  nonlinear  effects  of  the 
additional  lossy  modes  are  to  increase  the  output  coupling  and  to  both  broaden 
and  shift  the  laser  frequency. 

In  general,  the  solution  to  an  equation  with  a  nonlinear  source  term 
cannot  be  expressed  as  a  linear  superposition  of  solutions  to  the  homogeneous 
equation.  Only  when  one  of  the  resulting  expansion  coefficients  is  much  larger 
than  the  others  should  the  series  expansion  be  considered  as  a  numerically 
accurate  solution  to  the  nonlinear  problem.  If  one  is  interested  mostly  in 
quantitative  results,  the  method  of  using  series  expansions  in  empty  resonator 
modes  should  be  used  only  if  one  mode  dominates. 

However,  the  series  expansion  approach  can  always  be  used  to 
characterize  the  qualitatiave  behavior  of  high  energy  systems.  In  particular, 
this  approach  is  the  lowest  level  of  approximation  at  which  pulsations, 
bistability,  and  chaotic  behavior  involving  competition  between  transverse  modes 
can  be  expected  to  occur.  The  geometric  optics  approximation  cannot  lead  to  this 
type  of  unstable  behavior,  and  the  nonlinear  integral  equation  is  not  applicable 
in  the  region  of  chaos  or  pulsations.  Consequently,  the  equations  of  motion  of 
the  laser  using  series  expansions  in  empty  resonator  modes  may  be  the  best 
available  way  to  characterize  such  systems. 

We  have  not  considered  here  pulsed  mode  operation,  or  frequency 
dependent  gain.  Pulsed  operation  may  be  treated  by  allowing  the  expansion 
coefficients  to  have  a  slowly  varying  z  dependence,  which  adds  a  z  derivative  to 
the  equations  of  motion^^\  More  complicated  time  dependence  may  be  treated  by 
directly  solving  the  equations  of  motion  rather  than  simply  determining  the 
stable  points.  Frequency  dependent  gain  complicates  the  problem  because  the 
equations  are  no  longer  first  order  in  time  derivatives.  These  difficulties  have 
been  resolved  for  lasers  with  stable  cavities^^^^  and  therefore  the  same  methods 
should  be  applicable  for  unstable  resonator  cavities. 
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